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~~^This  report  presents  an  analytic  evaluation  of  the  differential 
Omega  experiment  being  conducted  by  the  Federal  Aviation  Administra¬ 
tion  in  the  Alaska/Yukon  region.  An  error  model  for  Omega  propaga¬ 
tion  prediction  residuals  i.s  used  to  predict  performance  achievable 
with  an  equipment  suite  consisting  of  a  TRACOR  599R  Omega  monitor 
receiver  and  associated  ground  station  equipment  built  by  Transport 
Canada,  and  an  airborne  suite  which  includes  a  modified  TRACOR  7620 
Omega  navigation  set  and  special-purpose  differential  Omega  equipment 
built  by  Systems  Control,  Inc.  (Vt.K  Predictions  are  made  for  use 
of  the  equipment  as  a  non-precision  approach  aid  within  the  differ¬ 
ential  Omega  coverage  area  and  for  use  as  an  enroute  navigation  aid 
between  areas  serviced  by  a  differential  Omega  ground  station. 
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INTRODUCTION 


1 . 


1.1  OVERVIEW 

The  Federal  Aviation  Administration  (FAA)  is  investi¬ 
gating  the  feasibility  of  using  differential  Omega  as  a  solu¬ 
tion  to  general  aviation  user  requirements  in  the  Alaska  region. 
At  present,  a  feasibility  demonstration  experiment  is  being 
conducted  in  the  Alaska/\’ukon  region.  The  experiment  is  a 
joint  venture  between  the  U.S.  Department  of  Transportation 
and  the  Canadian  Ministry  of  Transport. 

Differential  Omega  is  a  concept  for  improving  the 
navigation  performance  of  airborne  equipment  in  the  vicinity 
of  a  fixed-site  ground  station.  The  Omega  Navigation  System 
is  a  worldwide  Very-Low-Frequency  (VLF)  radio  navigation  aid 
(Ref.  1).  Because  the  errors  associated  with  Omega  are  highly 
correlated  in  space  (Ref.  2),  a  single  monitor  station  can 
service  a  relatively  large  area  (ranges  well  in  excess  of  100 
nm  from  the  monitor).  Basic  Omega  coverage  in  Alaska  is  good: 
signals  from  at  least  four  Omega  transmitting  stations  are 
normally  available  (Fef.  3).  Therefore,  differential  Omega 
offers  a  potential  cost-effective  alternative  to  establishing 
a  network  of  VOR/DME  stations  in  Alaska  (Refs.  3  and  4). 

The  current  differential  Omega  experiment  involves 
several  participants.  The  experiment  coordinator  is  the  Systems 
Research  and  Development  Serx'ice,  FAA.  Ground  station  equipment 


♦Very-High-Frequency  Omnidirectional  R.incr  (VOU)  .mil  Uist.incc- 
.Measuring  Equipment  (DME)  are  the  primary  aids  for  enroute 
navigation  in  the  continental  U.S. 
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was  built  by  Transport  Canada.  Airborne  equipment  is  installed 
in  FAA  aircraft.  Airborne  data  equipment  was  built  by  Sy.stenis 
Control,  Inc.  (Vt .  )  (SCI),  and  data  are  being  collected  by 
SCI.  The  airborne  navigation  receiver  equipment  was  manufac¬ 
tured  by  TRACOR ,  Inc. 


The  .Analytic  Sciences  Corporation  (TASC)  was  tasked 
by  the  FA.A  to  perform  an  analytic  evaluation  of  the  experi¬ 
ment  in  order  to: 


•  Determine  whether  the  experimental  differ¬ 
ential  Omega  system  can  meet  certification 
requirements  as  a  non-precision  approach 
aid 

•  Determine  the  suitability  of  the  Omega 
receiver  as  an  enroute  navigation  aid 

•  Define  primary  requirements  for  equipment 
suitable  for  certification  as  a  non¬ 
precision  approach  aid  and  for  enroute 
navigation . 


This  report  presents 


the  results  of 


these  analyses. 


1.2  DESCRIPTION  OF  THE  EXPERI.MENT 

The  differential  Omega  concept  is  illustrated  in  Fig.  1. 
Omega  navigation  signals  are  received  simultaneously  by  an 
airborne  receiver  and  a  ground  station  monitor.  At  the  ground 
station,  the  received  signal  is  compared  with  a  nominal  value 
and  a  correction  is  computed.  The  correction  is  uplinked  to 
the  aircraft  using  an  audio  signal  superimposed  on  a  Non-Direc- 
tional  Beacon  (NDB).  Airborne  equipment  decodes  the  telemetry 
message  and  applies  the  correction  to  improve  navigation  accuracy. 


selected 


For  the  Alaska/Yukon  experiment,  three  NDBs  have  been 
for  instrumentation  as  ground  stations:  CHENA  at 
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R.441&3 


Figure  1  The  Differenlial  Omega  Concept 

Fairbanks.  Alaska;  ISUX'IK  at  Inuv»ik,  Northwest  Territories; 
and  PUT  RIVER  at  Oeadhorse,  Alaska.  The  NDB  signal  has  a  usabl 
range  of  about  125  nm  and  the  relative  geometry  is  depicted  in 
Fig.  2 . 


Ground  station  equipment  consists  of  TRACOR  599R  Omega 
monitor  receivers,  a  specially  designed  signal  processor,  and 
recording  equipment.  The  receivers  have  been  modified  to  pro¬ 
vide  digital  output  of  phase  and  signal  amplitude,  and  atomic 
frequency  standards  are  used  as  reference  oscillators.  The 
differential  correction  message  format  is  summarized  in  Table  1 

Airborne  equipment  consists  of  an  interface  adapter/ 
message  decoder  for  receiving  the  correction  message,  a  TRACOR 
7620  Omega  navigation  set,  recording  equipment,  and  true  air¬ 
speed  and  stabilized  heading  inputs  to  the  navigation  set.  The 
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Figure  2  Geographical  Relalionship  of  the 
Ground  Stations 


TABLE  1 

DIFFERENTIAL  CORRECTION  MESSAGE  FOR>L^T 

'  '  ' 

Header 
Byte  Count 

Station  Identification 
Four  channels  of  10.2  kHz  data 
Four  channels  of  13.6  kHz  data 
Four  channels  of  11.33  kHz  data 
Signal  Quality  at  10.2  kHz 
Signal  Quality  at  13.6  kHz 
Signal  Quality  at  11.33  kHz 
Checksum 


axrborne  receiver  has  been  modified  to  accept  the  differential 

★ 

corrections  instead  of  computing  Omega  propagation  corrections  . 

1.3  REPORT  ORGANIZATION 

This  report  is  organized  into  five  chapters  and  three 
appendices.  Chapter  1  is  an  introduction  to  the  report  and 
provides  an  overview  of  the  experiment.  Chapter  2  presents 
the  essential  features  of  the  models  used  to  evaluate  the  per- 
of  the  TRACOR  7620  navigation  set  and  the  differential  Omega 
concept.  Details  of  the  Omega  models  are  described  in  Appen¬ 
dices  .A  and  B. 

.Analytical  predictions  for  the  differential  Omega 
system  are  presented  in  Chapter  3.  Section  3.2  describes  ap¬ 
plication  as  a  non-precision  approach  aid,  and  Section  3.3 
describes  application  as  an  enroute  aid.  Section  3.2  contains 
analyses  which  evaluate  the  sensitivity  of  the  performance 
predictions  to  uncertainties  in  model  parameters  and  modeling 
assumptions.  Additional  details  of  the  TRACOR  7620  mechaniza¬ 
tion  are  contained  in  Appendix  C. 

Chapter  ^  contains  an  analysis  of  the  performance 
predictions  in  terras  of  their  implications  for  certification 
of  the  differential  Omega  equipment.  The  results  are  present¬ 
ed  with  reference  to  three  separate  criteria:  certification 
requirements,  user  requirements,  and  procedural  requirements. 

The  criteria  are  defined  in  Section  4.1. 

Chapter  5  summarizes  the  report  and  presents  recom¬ 
mendations  for  further  work.  In  particular,  the  recommenda¬ 
tion  to  use  data  collected  at  the  ground  monitor  stations  to 
validate  the  models  of  Chapter  2  and  Appendices  A  and  B  is 
explained . 

*No  other  changes  have  been  made  to  the  TRACOR  7620  navigation 
algorithm,  and  TRACOR,  Inc.  did  not  have  the  opportunity  to 
optimize  the  design  for  the  differential  Omega  application. 
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ANALYTICAL  FORMULATION 


The  differential  Omega  experiment  will  involve  testing 
the  TRaCOK  7620  navigation  set  in  two  modes:  the  differential 
mode  when  monitor  corrections  are  available,  and  the  stand-alone 
.>mega  mode  when  corrections  are  not  available.  In  both  modes, 
the  navigation  set  must  perform  two  functions:  process  phase 
measurements  to  obtain  an  Omega  position  fix,  and  integrate 
additional  sensors  (viz.  true  airspeed  and  heading)  to  provide 
continuous  navigat ion  displays. 

Omega  position  fix  quality  is  evaluated  using  an  error 
model  for  Omega  signal  propagation  uncertainties  as  well  as 
for  receiver  errors.  The  Omega  propagation  model  described 
herein  is  a  statistical  model  derived  from  measurements  made 
on  the  Omega  Navigation  System  at  several  fixed  sites  through¬ 
out  the  world  (Ref.  2).  The  data  used  to  calibrate  the  model 
were  all  obtained  for  the  lowest  Omega  frequency,  10.2  kHz. 

.Application  of  this  model  to  this  differential  Omega 
experiment  relies  on  two  key  assumpt ions : 

•  The  same  model  structure  can  be  used  to 
describe  Omega  errors  at  13.6  kHz 

•  The  parameter  values  obtained  from  the 
worldwide  data  base  are  appropriate  for 
the  Alaska/Yukon  region. 

Parameter  values  obtained  from  Ref.  2  were  chosen  as  the  "nom¬ 
inal  values"  in  the  analysis,  and  an  analysis  was  performed  to 
evaluate  the  sensitivity  of  the  performance  predictions  to  vari¬ 
ations  in  the  parameter  values.  (The  sensitivity  analysis  is 


6 


presented  in  Chapter  3  in  the  form  of  curves  in  which  a  wide 
range  of  variation  about  the  nominal  parameter  value  is  evalu¬ 
ated.  ) 


Evaluation  of  the  TRACOK  7p20  as  a  navigation  set  in¬ 
volves  simulating  the  performance  of  the  navigation  algorithm 
along  specific  flight  profiles  in  the  Alaska/Yukon  region  using 
appropriate  environmental  parameters  for  that  region.  The 
simulations  are  performed  using  computer  programs  previously 
developed  by  TASC  for  the  TRACOK  7620  navigation  set  (Ref.  7). 
These  programs  model  the  TK,ACOR  7620  from  the  output  of  the 
receiver  phase  trackers  to  the  output  of  the  navigation  proc¬ 
essor. 


The  relationship  between  the  various  analytical  formu 
lalions  and  the  significance  attached  to  the  results  are  de¬ 
scribed  in  Section  2.1.  The  statistical  Omega  phase  error 
model  IS  summarized  in  Section  2.2.  (The  model  is  described 
in  more  detail  in  Appendix  A  and  a  theoretical  discussion  of 
the  most  significant  model  parameter,  the  correlation  coeffi¬ 
cient  between  errors  on  both  frequencies  from  the  same  Omega 
transmitting  station,  is  presented  in  Appendix  B.)  Section 
2.3  contains  a  brief  description  of  the  TRACOK  7620  navigation 
algor i t  hm . 


2 . 1  METHODOLOGY 

The  method  used  to  evaluate  the  performance  of  the 
differential  Omega  experiment  in  the  Alaska/^'ukon  region,  using 
the  TRACOR  7620,  is  summarized  in  Fig.  3.  Two  analysis  tech¬ 
niques  (statistical  and  direct  simulation)  are  used  to  analyze 
the  effects  of  three  error  sources:  Omega  phase  error,  wind 
error,  and  instrument  error.  The  contribution  of  each  of  these 
error  sources  to  navigation  error  in  the  TRACOR  7620  navigation 
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Figure  3  0\*erview  of  Analysis  Methodology 

filter  is  determined.  The  combined  contributions  of  the  indivi¬ 
dual  error  sources  provide  a  performance  estimate  of  the  accuracy 
of  the  TRACOR  7620  navigation  display.  In  addition,  accuracy 
obtainable  from  a  single  Omega  position  fix  is  evaluated  using 
two  techniques:  optimal  ( least  -  squares )  use  of  Omega  phase 
information  on  two  frequencies  defines  the  best  position  fixing 
available  to  a  navigation  set  with  a  similar  hardware  configura¬ 
tion,  and  use  of  the  specific  TRACOR  7620  composite  Omega  mechan¬ 
ization  describes  the  quality  of  the  input  to  the  navigation 
algorithm. 

Wind  errors  are  observ’ed  by  the  TRACOR  7620  navigation 
algorithm  through  the  Omega  position  fixes.  Changes  in  ground- 
speed  and  direction  result  from  shifts  in  the  wind,  since  the 
aircraft  normally  attempts  to  fly  at  constant  airspeed.  A 
direct  simulation  of  the  effect  of  winds  is  used  to  assess  the 
impact  of  wind  error.  The  computer  programs  used  in  the  analy¬ 
sis  were  previously  developed  by  TASC  for  the  TRACOR  7620  navi- 
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gation  set  (Kef.  7).  Satisfactory  filter  response  to  a  worst- 
case  wind  error  condition  is  used  to  establish  the  suitability 
of  the  filter  for  this  application.  Additional  information  is 
supplied  to  the  navigation  filter  by  a  heading  reference  and 
an  airspeed  indicator.  Reliable  error  models  for  these  instru¬ 
ments  were  not  available  for  this  study,  and  the  contribution 
of  nominal  errors  from  these  instruments  to  navigation  error 
is  not  considered.  However,  the  impact  of  instrument  fai lure 
on  navigation  error  is  studied  by  direct  simulation  techniques. 
The  combined  effect  of  Omega  phase  error,  wind  error,  and  in¬ 
strument  failure  on  the  TRACOR  7b20  navigation  filter  provides 
an  estimate  of  the  accuracy  of  the  total  navigation  display  in 
the  aircraft.  Since  nominal  instrument  errors  are  not  con¬ 
sidered,  accuracy  predictions  for  the  total  navigation  display 
under  nominal  conditions  apply  only  to  a  navigator  with  perfect 
heading  and  airspeed  information , 

To  evaluate  the  full  potential  of  the  differential 
Omega  concept,  an  analysis  is  performed  to  determine  the  accu¬ 
racy  achievable  from  usage  of  phase  information  in  a  single 
Omega  position  fix.  A  lower  bound  on  Omega  position  fix  accu¬ 
racy  is  provided  by  opt imal  ( least  -  squares  )  processing  of  the 
dual  -  frequency  phase  information  from  four  stations.  A  second 
measure  of  Omega  position  fix  accuracy  is  provided  by  the  com¬ 
posite  mechanization  used  as  an  input  to  the  TRACOR  7620  navi¬ 
gation  filter.  The  composite  mechanization  evaluation  does 
not  directly  measure  the  performance  of  any  part  of  the  TRACOR 
7620  filter:  however,  it  does  provide  a  lower  bound  on  the 
Omega  position  fix  accuracy  used  as  an  input  to  the  filter. 

The  performance  of  these  two  estimators  (optimal  and  composite) 
serves  to  establish  the  performance  obtainable  with  differential 
Omega  operating  as  a  simple  disp lay  of  instantaneous  position 
without  considering  the  navigation  function.  Although  the 
potential  benefits  associated  with  phase  measurement  averaging 
are  not  realized,  the  errors  and  response  characteristics  asso- 
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cidlrd  With  a  specific  navigation  mechanization  are  also  elimi* 
nated.  Omega  position  fix  evaluation  also  provides  a  convenient 
mechanism  to  evaluate  the  sensitivity  of  performance  to  changes 
in  parameters  of  the  Omega  phase  error  model. 

2.2  PH.ASE  ERROR  .MODEL 

Navigation  error.  is  defined  as  the  vector  differ¬ 

ence  between  the  rece i ver - i nd i ca t ed  position,  x^ .  and  the  true 
pooitiori,  x^  ,  where: 


^  =  X j  -  x^  ( 1  ) 

It  IS  convenient  to  express  navigation  error  in  terms  of  north 
and  east  position  error  components  on  the  surface  of  the  earth 
(dN  and  .  respectively).  If  the  position  deviations  are 
small  relative  to  the  range  to  an  Omega  transmitting  station, 
the  equivalent  phase  dev'iation  associated  with  the  position 
deviation  is  a  simple  function  of  the  bearing  angle  to  the 
station.  Omega  phase  measurements  also  include  errors  caused 
by  propagation  anomalies  and  receiver  measurement  errors.  The 
resulting  linearized  scalar  expression  for  the  measured  phase 
deviation  caused  by  navigation  error  and  by  propagation  and 
receiver  errors  is 


cos  t  ^ 


AN  •» 


sin  t.  AE  ■* 


where 


(2) 


is  the  measured  phase  deviation  from 
station  i  on  frequency  j 


\.  is  the  wavelength  of  the  signal  with 
frequency  j 

is  the  bearing  angle  to  station  i 
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6o.  .  is  the  phase  error  t rom  station  i  on 
^  frequency  j  caused  by  propagation 
error  and  receiver  error. 

It  is  convenient  to  write  Eq.  2  in  matrix  form  by  defining  a 
vector  of  phase  measurement  deviations,  which  includes  as 

components  all  of  the  eight  available  signals: 

^  =  H  ^  ^  (3) 

where 

H  is  the  measurement  geometry  matrix 

is  the  vector  of  north  and  east  posi- 
t ion  errors 

is  the  vector  of  phase  errors. 

When  signal  phase  measurements  from  more  than  three 
separate  transmitting  stations  are  made,  redundant  information 
is  available  which  can  be  used  to  reduce  navigation  error. 

The  impact  of  Omega  phase  error  on  navigation  error  depends  on 
the  manner  in  which  the  navigation  set  uses  the  phase  measure¬ 
ments.  As  previously  noted,  two  techniques  for  using  phase 
information  are  evaluated: 

•  Optimal  -  A  minimum-variance  (*'least- 
squares")  estimate  of  position  based  on 
all  available  phase  measurements 

•  Composite  Mechanization  -  A  suboptiroal 
estimate  of  position  based  on  the  speci¬ 
fic  combination  of  phase  measurements 
implemented  at  the  input  to  the  TRACOR 
7620. 

Note  that  the  optimal  formulation  will  always  provide  the  smal¬ 
ler  error  (under  nominal  conditions):  a  suboptimal  technique 
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such  as  the  TRACOR  7620  mechanization,  however,  may  be  simpler 
to  implement  and  may  possibly  exhibit  less  sensiliv'ity  to  mod¬ 
eling;  errors. 

2.2.1  Optimal  Processing  of  Phase  Measurements 

The  minimum  variance  estimate,  of  the  position 

deviations  is  (Ref.  10> 


(^) 


where 

R  =  El^ 

E|-l  denotes  the  expected  value  operator. 


The  covariance  of  the  error  in  the  estimate  ^  is 

P  =  (H^R'^H)’^  (5) 

where 

P  =  E{ 

A  statistical  model  for  is  required  to  form  R.  Details 

of  the  model  for  «re  presented  in  Appendix  A.  A  summary 

of  the  salient  characteristics  of  the  model  is  presented  in 
Table  2. 

2.2.2  Composite  Mechanization 

In  the  TRACOR  7620  navigation  set,  eight  Omega  phases 
are  combined  to  form  three  composite  phase  measurements  with 
the  following  general  form: 
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TABLE  2 

SUMMARY  OF  COMPONENTS  OF  NOMINAL  PHASE  ERROR 
covariance  matrix  MODEL 


Z-iUit 


kASljOM  HROHACATlOk  F.KKOh  NOOtL 
‘•''■*1  2‘-l  -  ‘l  ** 


*  CAP  !  .  !  -  I  Cl  ,*  r»p  J  — -  -  ii  .  o/  • 


.  o/  .  |l2  -  .j|j 

*  l‘2  -  'll)j 

li  =  1 SOO  iia  Oj*  =  9  83  CfK^  j  -  ijj  =  jvrrijje  •untlor  upd<<(«  deljy  -  3  Mn 


I  j  =  0  9^.3  hi 

r.  s  12  hi 


:  8  hr 


u  ,*  t  17  8  Cft* 

»  h  13  cec^ 

t  ‘9 

O  •  *  I  4b  cec* 


EIA*’’ 


1.1*  b*^2  • '  2 I  .  1  3  h*^P*pl  *  '■’*  *•'**!  .10. 2*-2’'2*  **!,10  2'-l''l 


E|/.»  ,1  j  )  I  s  |,w  10  2*52'*2***i  .10  2'-l’‘l 


( X  j  .  t  j  )  I  .  W  s  1 

(», .1, ) I .  p  »  0  5 


tclCLIVf.R  tKROh  flODr.l 


r.iA*  ^  ^ 


*ij.c  ^  ■  *j.l3.6^  ‘  ^^®i.l0.2  ‘  ®j,10.2^ 


where 


•i ,13.6 

*j.l0.2 

a,b 


is  the  composite  phase  from  stations 
i  and  j 

is  the  phase  from  station  i  at  13.6  kHz 

is  the  phase  from  station  j  at  10.2  kHz 

are  constants  (a  =  lS/8  and  b  ~  3/2  for  the 
TRACOR  7620). 
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The  relation  between  composite  phase  deviation  and  position 
error  for  this  set  is  (from  Eqs.  2  and  f>) 


'ij.c  '  * 


^(sin  f  .-sin 


ij  .c 


where  X  is  the  composite  wavelength  given  by  (r — - —  z — - — \ 

^  \^13.6  ^10.2/ 


If  the  scalar  composite  phase  deviations  are  combined  to  form 
a  vector  of  composite  phases  (denoted  by  ,  the  derivation 
of  Eqs.  3-5  can  be  applied  to  the  TRACOR  7620  processing. 


2.3  TRACOR  7620  NAVIGATION  FILTER 


The  TRACOR  implementation  of  the  Omega  navigation 
algorithm  is  a  suboptimal  Kalman-like  filter  (Ref.  17).  It 
contains  a  four-state  dynamic  model  that  allows  position  pre¬ 
diction  and  correction.  This  section  presents  a  derivation  of 
an  error  model  for  this  algorithm  and  describes  some  of  its 
characteristics . 


TRACOR 


lorithm 


The  TRACOR  algorithm  has  two  sections;  a  dead- reckoning 
section,  which  uses  airspeed  and  heading  information,  and  an  up¬ 
date  section  that  processes  Omega  measurements. 

Dead-Reckoning  -  The  dead- reckoning  section  performs 
a  rectangular  integration  of  a  four-element  state  vector  every 
1.25  sec.  The  filter  states  are  latitude  and  longitude,  in 
radians,  and  latitude  and  longitude  wind  rate,  in  radians  per 
second : 
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(8) 


X  =  ILAT,  LON.  WINDj^. .  VlNDj.)^ 
Tht-  dead- reckoning  equation  is: 


x(kT)  = 


'l 

0 

1.25 

0 

cos  TH 

0 

1 

0 

1.25 

x(kT-T>-^l  .25 

sin  TH 
cos  LAT 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

TAS 

3600R 


<9) 


where 


TH  IS  the  true  heading 
TAS  IS  the  true  airspeed  (kt) 


K  IS  the  radius  of  the  earth  (nm) 
e 


Omega  Updates  -  Every  10  sec.  the  algorithm  processes 
a  set  of  Omega  measurements  and  applies  a  correction  to  the 
current  state  estimate: 


xikl")  =  x{kT')  ♦  K|z  -  «^(x(kT‘))l  (10) 

where  z  is  the  vector  of  observ'ed  phases  and  ♦^(x)  is  the  vec¬ 
tor  of  estimated  composite  phases,  as  a  function  of  estimated 
pos 1 t ion . 


The  filter  gain  matrix  K  is  variable  and  depends  on 
estimated  position: 


K  =  Pj  H(x)  (11) 

and  are  fixed,  while  H  depends  on  station  geometry. 

The  TRACOR  7f>20  navigation  filter  does  not  update  the  matrix 
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2.3.2  TRACOR  Error  .Model 

From  Appendix  C,  the  position  and  wind  errors  in  the 
TRACOR  7620  state-vector  estimates  are  given  by 

^(kT)  =  11-KHJ  !♦  ^(kT-T)  ♦  L(w(kT-T)  -w(kT)) 

GV|j^(kT)l  ♦  Kv(kT)  (12) 

where 


$  IS  the  state  transition  matrix 

C  is  the  noise  weighting  matrix 

L  IS  the  control  weighting  matrix 

w  IS  the  wind  vector 

V  IS  the  instrument  error  vector 

-m 

V  IS  the  Omega  phase  measurement  error  vector. 


Equation  12  is  used  to  evaluate  the  steady-state  response  of 
the  TRACOR  7620  navigation  filter  to  these  error  sources. 


3. 


perform/\nce:  predictions 


3.1  INTRODUCTION 

Differential  Omega  is  under  consideration  as  a  non- 
precision  approach  aid  and  as  an  enroute  aid  between  airports. 

In  this  chapter,  the  model  for  differential  Omega  phase  error 
presented  in  Appendix  A  and  the  TRACOR  7620  navigation  set 
implement  at  ion  described  in  Appendix  C  are  used  to  predict 
navigation  performance  for  both  applications.  In  Section  3.2, 
the  expected  nominal  position  error  near  an  airport  is  pre¬ 
sented  along  with  the  performance  sensitivity  to  assumed  operat¬ 
ing  conditions  and  to  the  .assumed  model  parameters.  The  use 
of  differential  corrections  in  the  enroute  mode  is  analyzed  in 
Sect  ion  3.3. 

3.2  USE  AS  A  NON-PRECISION  APPROACH  AID 

3.2.1  Nominal  Performance  < 


The  expected  performance  of  a  differential  Omega  navi 
gation  set  on  an  approach  to  Inuvik,  2  nm  away,  is  presented 
in  Table  3.  Because  of  the  high  sensitivity  of  the  composite 
mechanizat ion  to  the  correlation  between  frequencies,  p,  (see 
Section  3.2.3)  results  are  presented  for  a  worst-case  nominal 
correlation,  p=0.5,  and  also  for  p=0.9.  In  addition,  perform¬ 
ance  predictions  for  a  single- frequency  mechanization  as  well 
as  performance  predictions  with  no  differential  corrections 
are  included  for  comparison.  Along-  and  cross-track  position 
error  estimates  are  provided  for  an  approach  to  Inuvik  at  a 
geodetic  bearing  of  90  deg.  The  along-  and  cross-track  esti¬ 
mates  in  Table  3  are  combined  to  give  radia’  position  error  in 
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Table  ‘4.  Kadial  position  error  is  used  as  the  measure  of  per 
formance  throughout  the  remainder  of  this  chapter. 

Results  of  the  position-fix  error  analysis  show: 


•  Optimal  posit  ion- fix  performance  is  not 
sensitive  to  the  correlation  coefficient 
p.  The  total  error  of  0.15  nm  is  produced 
almost  entirely  by  propagation  prediction 
errors 

•  Nominal  composite  mechanirat ion  position- 
fix  perform.ince  for  0  =  0.*^  has  a  position 
error  due  to  propagation  errors  which  is 
nearly  optimal,  but  a  relatively  large 
contribution  of  receiver  error  causes 
performance  to  be  significantly  inferior 
to  the  optimal  pos  1 1  ion- f  i.\  result 

•  f’rocessing  phase  measurements  from  only 

1 3 . b  kHz  signals  provides  neatly  optimal 
posit  ion- fix  performance  which  is  signi¬ 
ficantly  superior  to  the  compos i t e - f re - 
quency  mechanization 

•  Differential  corrections  from  a  monitor 
are  effective  in  reducing  position  error 
in  all  cases. 


3.2.2  Sensitivity  to  Operating  Conditions 


In  thi.s  section,  the  impact  on  navigation  accuracy  of 
changes  in  the  following  operating  conditions  is  presented: 


•  .'\irport  location 

•  Vehicle  maneuvers 


Station  outage 
Loss  of  instruments 

Monitor  update  delay  (signal  averaging 


Sensitivity  lo  Airpori  Location  -  Since  the  estima¬ 
tion  gain  matrix  depends  on  station  receiver  geometry, 
position  error  response  uill  change  from  one  location  to  another 
In  particular,  response  improves  as  the  vectors  to  the  stations 
become  more  nearly  perpendicular  (Ref.  7).  Nevertheless,  sta¬ 
tion  geometry  does  not  vary  greatly  in  the  region  of  interest. 
(See  Fig.  4.  Arrows  represent  range  and  bearing  to  transmit¬ 
ting  stations.  .M.M  denotes  megameter.)  Geometry  appears  to  be 
the  best  at  Fairbanks  and  the  worst  at  Inuvik. 

Table  5  shows  the  effect  of  geometry  on  the  statis¬ 
tical  performance  of  the  position-fix  and  navigation  filter 
position  estimates.  Performance  at  locations  near  the  three 
different  airports  under  consideration  is  nearly  the  same,  but 
IS  generally  best  at  Fairbanks  and  worst  at  Inuvik. 

To  test  the  effect  of  these  different  geometries  on 
the  TKACOK  7620  navigation  filter  dynamics .  simulated  error 
histones  were  generated  for  areas  near  Fairbanks  and  Inuvik. 


TABLE  3.2-3 

R.MS  R.AI’IAL  position  ERROR  FOR  DIFFERE.N'T  AIRPORTS 
(2  nm  FROM  AIRF’ORT) 


RMS  RADIAL  POSITION 

ERROR  (nm) 

1 KPORT 

INTAIK 

FAIRBANKS 

DEADHORSE 

CATEGORY 

..  .   .i 

Posi t ion- f  ix 

Composite  mechanization 
(p=0.5) 

0.40 

0.35 

0.36 

Optimal  (p=0.5) 

0.15 

0.13 

0.14 

TR.-\C0R  7620 

Navigation  Filter 


0.33 


0.34 


0.34 


«« t*>  • 


<1  ni*t»HO*«si 


bl  IMUVI*. 


cl  r*lM»AMKS 


Figure  i*  Omega  Geometry  at  Selected  Sites 


Only  three  Omega  signals  were  simulated  in  order  to  accentuate 
the  effects  of  geometry.  All  other  conditions  were  held  con¬ 
stant  : 

•  Airpseed:  250  kt 

•  Initial  track:  180“ 

•  Wind  steps  from  zero  to  10  kt  after  2 

min  into  flight 

•  Aircraft  turns  clockwise  through  270“ 

(from  0“)  at  3“/sec  after  10  min  into 
flight 

•  Signals  from  transmitting  stations 
C,  D,  H  assumed  to  be  available 
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•  Noise-free  reception  and  trror-free  dif¬ 
ferential  corrections  available 

•  Continuous  true  airspeed  and  stabilized 
heading  data  available. 

Figure  5  depicts  the  simulation  results  for  trajectories  near 
Fairbanks  and  Inuvik.  Maximum  position  error  for  the  Inuvik 
simulation  is  approximately  10»  greater  than  the  maximum  posi¬ 
tion  error  for  the  Fairbanks  simulation.  Error  dynamics  are 
essentially  the  same  at  both  locations. 

These  results  demonstrate  that  Omega  geometry  has 
very  little  effect  on  pe rformance  over  the  region  of  interest . 
However,  performance  near  Inuvik  may  be  taken  as  worst-case. 

Maneuver  Sensitivity  -  Since  the  TRACOK  navigation 
filter  depends  on  a  1.25  sec  rectangular  integrator,  it  may  be 
possible  for  a  maneuver  to  create  unacceptably  large  position 
errors.  A  simulated  error  history  was  generated  to  determine 
the  sensitivity  of  position  error  to  a  maximum  procedure  turn 
for  the  following  conditions; 

•  Initial  position  hH“S.  13<«®W 

•  Airspeed:  250  kt 

•  Initial  track.  180* 

•  Aircraft  turns  clockwise  through  270*  at 

3* /sec  after  80  min  into  ^lignt 

•  Signals  from  t  r ‘.nsm  1 1 1  mg  stations 
A.  C.  D,  H  assumed  to  be  available 

•  Noise-free  reception  and  error-free  dif¬ 
ferential  corrections  available 

•  Continuous  true  airspeed  and  stabilized 
heading  data  available. 

Results  of  this  simulation  are  shown  in  Fig.  6. 
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Figure  5  Effect  of  Airport  Location  on  TRACOR  7620 
Navigation  Performance 
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Effect  of  Maneuvers  on  TRACOR  7620 
Navigation  Performance 


Initially,  small  position  errors  occur  due  to  numer* 
ical  errors  in  simulation  initialization.  At  the  start  of  the 
turn,  a  sharp  transient  appears.  It  reaches  a  maximum  radial 
position  error  of  approximately  0.02  nm  and  quickly  decays. 
Athough  a  maneuver  introduces  some  position  error,  the  resul¬ 
ting  maximum  radial  error  is  not  likely  to  exceed  0.02  nm. 

The  maneuver- induced  error  decays  quickly,  and  system  stabil¬ 
ity  is  not  compromised  by  maneuvers. 


Sensitivity  to  Station  Outage  -  In  the  Alaska/Yukon 
region,  signals  for  Omega  transmitting  stations  A,  C,  D,  and  H 
can  generally  be  received.  In  the  event  of  signal  loss  from 
one  of  these  stations,  performance  degrades  to  the  levels  shown 
in  Table  6.  Performance  degradation  is  largest  for  loss  of 
signals  from  either  station  A  or  H. 


TABLE  b 

RMS  RADIAL  POSITION  ERROR  FOR  STATION  OITAGE  CONDITIONS 

(2  nm  FROM  1NU\’1K) 


To  examine  the  effect  of  the  loss  of  one  station  on 
the  position  error  dynamics  of  the  TRACOR  7620  filter,  two 
simulations  were  run.  Figures  5a  and  5b  present  results  for  a 
receiver  near  Inuvik  when  signals  from  transmitting  stations 
C,  D,  and  H  are  received.  In  Fig.  7  results  correspond  to 
signals  from  stations  A,  C,  D  and  H  being  used.  Other  Simula* 
tioii  parameters  are: 


•  Initial  position:  68®N,  134®V 

•  Airspeed:  250  kt 

•  Initial  track:  180® 

•  North  wind  steps  from  zero  to  10  kt  after 
2  min  into  flight 

•  Aircraft  turns  clockwise  through  270®  at 
3®/sec  after  20  min  into  flight 

•  Noise-free  reception  and  perfect  differ¬ 
ential  corrections  available 

•  Continuous  true  airspeed  and  stabilized 

heading  data  available. 


26 


*  «4I^I 
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d)  North  Position  Error  b)  East  Position  Error 

Figure  7  Effect  of  Signal  Availability  on  TRACOR  7620 
Navigation  Performance 

The  maximum  radial  position  error  increases  50%  with  this  loss 
of  station  A,  from  approximately  0.2  nm  with  signals  from  four 
stations  received,  to  0.3  nm  with  signals  from  three  stations 
received.  Also  note  that  the  maximum  radial  error  occurs  at  a 
slightly  later  time  in  the  three-station  case  than  in  the  four- 
station  case. 

Sensitivity  to  Loss  of  Instruments  -  The  TRACOR  7620 
position  dead- reckoner  requires  airspeed  and  heading  informa¬ 
tion.  If  these  inputs  are  lost,  position  errors  will  increase. 
Loss  of  heading  and  airspeed  information  are  simulated  separ¬ 
ately.  as  described  below. 

First,  assume  that  true  heading  data  are  lost  and  the 
filter  continues  to  use  the  last  known  heading.  As  long  as 
the  aircraft  continues  to  fly  in  a  straight  line,  no  errors 
occur.  However,  large  transient  position  errors  arc  expected 
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during  turns.  Eventually,  the  filter  will  estimate  apparent 
wind  to  compensate  for  the  fact  that  the  aircraft  is  flying  in 
the  "wrong"  direction. 

To  test  the  effects  of  loss  of  heading  information,  a 
simulation  was  run  with  the  following  parameters: 

•  Initial  position:  bS°N, 

•  .■\irspeed:  250  kt 

•  Initial  t  rack :  1 80“ 

•  Constant  heading  of  180“  applied  to 
filler 

•  .Aircraft  turns  clockwise  through  270“  at 

3“ /sec  after  10  min  into  flight 

•  Signal*  from  transmitting  stations  A,  C. 

D.  and  H  assumed  to  be  available 

•  Noise-free  reception  and  error-free  mon¬ 
itor  corrections  assumed 

•  Continuous  true  airspeed  data  available. 

.^s  Fig.  8  shows,  large  transient  errors  occur  and  reach  a  ma.x- 
imu.m  radial  position  error  of  approximately  8.5  nm.  This  de¬ 
terministic  effect  is  sufficiently  large  to  dominate  the  random 
errors  associated  with  Omega  phase  measurements.  Recovery  of 
TR.^COR  7620  filter  takes  more  than  20  min  since  measurements 
must  propagate  through  two  slate  variables,  wind  and  position. 
.^Iso  note  that  north  and  east  wind  estimates  reach  steady-state 
values  of  250  kt .  This  would  be  sufficient  wind  to  move  the 
aircraft  east  at  the  known  true  airspeed  even  if  the  true  heading 
were  south,  as  the  dead- reckoner  assumes. 

Now,  suppose  that  airspeed  information  is  lost  and 
the  filter  continues  to  use  the  last  known  value,  2-transform 
analysis  (see  Appendix  C)  indicates  that  long-term  position 
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Figure  8  Effect  of  Loss  of  Heading  Data  on 
TRACOR  7620  Navigation  Performance 


errors  due  to  increasing  airspeed  errors  should  be  slowly- 
growing  but  unbounded.  Wind  estimates  should  compensate  for 
most  of  the  airspeed  error.  To  test  the  effect  of  loss  of 
airspeed  data  on  TRACOR  7620  position  error,  a  simulation  was 
run  with  these  parameters: 


29 


•  Initial  position:  b8“N ,  13^°V 

•  Initial  airspeed:  250  kt 

•  Airspeed  decreases  continuously  by  100 

kt/hr 

•  Constant  airspeed  data  applied  to  filter 

•  Aircraft  turns  clockwise  through  270“  at 
3° ^sec  after  10  min  into  flight 

•  Signals  from  t ransna  1 1 i ng  stations  A.  C, 
D,  and  H  assumed  to  be  available 

•  Noise-free  reception  and  perfect  differ¬ 
ential  corrections  assumed 

•  Continuous  stabilized  heading  data 
avai 1 abl e . 


Figure  9  shows  the  s 1 owly - growi ng  position  errors.  Note  that 
the  position  error  transient  starts  at  10  min,  as  the  aircraft 
turns.  The  apparent  headwind  is  forced  to  shift  as  the  heading 
changes.  Also  note  that  the  estimated  wind  ramp  has  a  slope 
of  approximately  100  kt/hr,  as  the  filter  compensates  for  in¬ 
creasing  indicated  airspeed  error. 


The  result  of  loss  of  instruments  can  be  sumitiarized 
as  follows: 


•  Loss  of  stabilized  heading  data  results 
in  large  position  error  transients  after 
each  maneuver  that  eventually  damp  to 
zero 

•  Loss  f  airspeed  information  results  in 
relatively  slowly-growing  position  errors 
which  eventually  become  unbounded. 
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Effect  of  Loss  of  Airspeed  Data  on 
TRACOR  7620  Navigation  Performance 


Update  delay  -  Under  nominal  conditions,  phase  meas 

★ 

urements  at  the  monitor  are  averaged  for  6  min  before  trans 


♦Averaging  is  a  consequence  of  the  narrow  tracking  bandwidth 
in  the  TRACOR  599R  monitor  receiver.  The  filters  being  used 
by  Transport  Canada  for  this  experiment  have  an  equivalent 
averaging  time  of  approximately  5  min. 
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mission  to  the  navigator  to  rt-duce  t  ht*  t-lfecl  of  noise  at  the 
monitor  receiver.  Since  the  random  component  of  propagation 
error  changes  with  time,  an  error  is  i nt  roduced  bv  averaging 
because  of  the  resulting  "staleness"  of  the  monitor  update 
with  respect  to  the  phase  measurement  at  the  navigator.  Fig¬ 
ure  10  shows  that  as  the  averaging  time  is  increased  from  zero 
position  error  first  decreases  because  of  monitor  noise  aver¬ 
aging.  but  then  increases  due  to  staleness  of  the  update.  The 
optimal  averaging  time  appears  to  be  about  1  min.  Note  that 
tne  effective  update  delay  (ltT-tj|in  the  model.  Table  2),  is 
one-half  of  the  signal  averaging  time. 

3.2.3  Sensitivity  to  .Model  Parameters 


The  statistical  model  used  to  assess  the  performance 
of  the  differential  Omega  concept  in  the  Al aska/\’ukon  region 
IS  based  on  limited  data  collected  worldwide.  Since  the  model 
parameters  appropriate  for  the  .Alaska/Yukon  region  may  be  some 
what  different  than  those  in  the  worldwide  model,  it  is  impor¬ 
tant  to  assess  the  impact  of  changes  in  model  parameters  on 
navigation  accuracy.  In  this  section  the  effect  of  changes  in 
the  following  model  parameters  is  considered: 

•  Correlation  coefficient  between 
f  requenc i es  ( p  ) 

•  Correlation  distance  (d) 

•  Correlation  time  ( t j  ) 

•  Standard  deviation  of  random  phase  error 

(/pj"  4  02^  ♦  ♦  0^^) 

•  Receiver  error. 

Correlation  coefficient  -  As  shown  in  Fig.  11,  an 
increase  in  the  frequency  cross-correlation  coefficient  above 


32 


I  A 


X  I 

4 


Sens 
to  S 
( 2  n 


the  (conservative)  nominal  of  p=0.b  indicates  improved  perform¬ 
ance  for  the  composite  implementation  position-fix  and  TRACOR 
7620  navigation  filter.  Notice,  however,  that  even  when  p=l , 
performance  using  only  13.6  kHz  signals  is  superior  to  the 
composi te- f requency  implementation  due  to  the  contribution  of 
receiver  error.  (See  Table  U.) 

Correlation  distance  -  Figure  12  shows  that  sensi¬ 
tivity  to  the  correlation  distance  parameter,  d,  in  the  model 

Ts  ~?.\T'feii5riy  sm'ai“i  njT~t?’‘jOti“rtnfT — T+re  sriis'i'c  . — 

rapidly  as  the  assumed  correlation  distance  approaches  the 
monitor-navigator  separation  of  2  nm. 

Correlation  time  -  Figure  13  shows  the  sensitivity  of 
navigation  performance  to  the  correlation  time,  tj.  Perform¬ 
ance  IS  very  sensitive  to  a  decrease  in  tj  as  it  approaches 
the  nominal  update  delay  time  of  3  min. 

Standard  deviation  of  propagat ion- i nduced  phase  error  - 
If  the  receiver  error  is  assumed  to  be  zero,  then  rms  radial 
position  error  is  a  linear  function  of  the  standard  deviation 
of  propagat ion- 1 nduced  phase  error.  Figure  1^  shows  this  linear 
function  combined  with  the  nominal  receiver  error  contribution 
to  position  error  to  produce  total  rms  radial  position  error 
as  a  function  of  propagat ion- induced  phase  error  standard  devi¬ 
ation.  The  sensitivity  increases  as  the  residual  propagation 
error  increases. 

Receiver  error  -  If  the  roles  of  propagat ion- induced 
error  and  receiver  error  are  reversed  in  Fig.  14,  the  sensitiv¬ 
ity  of  rms  radial  position  error  to  receiver  error  is  obtained. 
Results  are  shown  in  Fig.  15.  As  expected,  the  position  error 
increases  with  an  increase  in  receiver  error. 
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Figure  15  Sensitivity  of  RMS  Radial  Position  Error 
to  Receiver  Error  (2  nm  from  Inuvik) 
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3.3 


USE  AS  AN  ENROUTE  AID 


In  the  previous  section,  performance  of  the  differen¬ 
tial  Omega  system  on  an  approach  to  Inuvik  was  evaluated.  In 
this  section,  the  performance  of  differential  Omega  is  evalu¬ 
ated  for  a  flight  away  from  Inuvik  (and  consequently,  away 
1 rom  a  monitor  site).  Section  3.3.1  describes  navigation  per¬ 
formance  while  the  aircraft  is  within  range  of  the  monitor 
update  signal.  Section  3.3.2  evaluates  the  alternatives  when 
a  differential  update  is  no  longer  receivable.  Section  3.3.3 
provides  a  summary  of  recommended  enroute  uses  of  differential 
Omega . 


It  is  assumed  that  a  differential  correction  can  be 
received  up  to  125  nm  from  the  monitor.  Figure  16  shows  navi¬ 
gation  performance  as  a  function  of  distance  from  the  monitor, 
up  to  125  nm.  Performance  degrades  because  of  spatial  decor¬ 
relation  of  the  monitor  information,  but  the  differential  cor¬ 
rection  still  significantly  improves  performance. 

3.3.2  Enroute  Performance  Beyond  Monitor  Range 

When  a  differential  update  can  no  longer  be  received, 
three  operational  options  are  considered: 

Option  1:  Apply  the  last  differential  update  re¬ 
ceived  as  if  it  were  current. 

Option  2:  Use  the  last  differential  correction 
to  calibrate  the  onboard  Predicted 
Propagation  Corrections  (PPCs)  at  the 
monitor  site  and  use  the  calibration 
to  adjust  PPCs  at  the  indicated  air¬ 
craft  position. 

Option  3:  Ignore  the  monitor  information  and 
navigate  using  PPCs. 
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Figure  Jb  Savigiition  Performance  Within  Range 
of  Monitor  Updates 


Performance  under  these  three  options  is  evaluated  below: 

Option  1  -  There  are  attractive  reasons  for  continu¬ 
ing  to  apply  the  last  received  differential  correction.  Such 
a  policy  IS  easy  to  implement,  since  it  does  not  require  PPC 
calculations.  PPCs  are  knouTi  to  be  highly  correlated  in  space , 
suggesting  that  the  differential  correction  also  contains  use¬ 
ful  information  for  some  distance  from  the  monitor  recei\»er. 
However,  PPCs  exhibit  rather  significant  temporal  variations. 

To  test  the  effect  of  temporal  PPC  variations  on  position  error, 
a  simulation  was  run  with  the  following  conditions: 


•  Initial  position:  64.5®N,  1A7.5®V 

•  Airspeed:  250  kt 
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Track:  0*' 


•  Signals  from  transmitting  stations  A,  C, 

D>  and  H  assumed  to  be  available 

•  Noise- free  reception  assumed 

•  Include  phase  shifts  corresponding  to 
2100-2200  GMT.  23  December 

•  Assume  that  the  last  differential  correc¬ 
tion  is  obtained  at  the  initial  time  and 
position 

•  Assume  that  the  filter  does  not  apply 
PPCs. 

The  simulated  conditions  characterize  very  large, 
probably  worst  case,  PPC  changes.  There  are  other  dates  and 
times  when  the  simulated  effect  is  minor.  However,  it  is 
clear  from  Fig.  17  that  this  option  for  area  navigation  can 
give  rise  to  unacceptably  large  errors,  regardless  of  other 
effects  present. 


TiMf  Immi  TIME  (mml 

a)  North  Position  Error  b)  East  Position  Error 

Figure  17  Effect  of  Uncompensated  PPC 
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Option  2  -  This  option  requires  the  Omega  navigation 
Set  to  compute  two  PPCs  for  each  differential  Omega  measure¬ 
ment:  one  for  the  estimated  current  position  of  the  aircraft, 

and  one  for  the  location  of  the  monitor.  By  removing  the  mon¬ 
itor  PPC  from  the  differential  update  and  replacing  it  with 
the  aircraft  PPC.  performance  degradation  caused  by  changes  in 
the  PPCs  between  the  monitor  and  the  aircraft  is  avoided. 


Figure  lb  (solid  curves)  shows  that  this  technique  provides 
improved  performance  for  up  to  1.5  hr  after  the  last  differen¬ 
tial  update  compared  to  performance  with  no  differential  up¬ 
date.  This  improvement  is  possible  because  of  the  relatively 
long  temporal  correlation  of  residual  propagation  error  informa 
tion  in  the  monitor  update.  Assuming  an  aircraft  speed  of  250 
kt ,  Fig.  19  (solid  curves)  shows  that  monitor  information  pro¬ 
vides  improved  performance  up  to  500  nm  from  the  monitor.  The 
dashed  curve  in  Fig.  18  shows  performance  in  the  hypothetical 
case  that  temporal  decorrelation  alone  is  present,  which  is 
equivalent  to  assuming  that  the  aircraft  remains  125  nm  from 
the  monitor.  In  Fig.  19  the  dashed  curve  shows  performance 


assuming  that  only  spatial  decorrelation  occurs,  which  is 
equivalent  to  assuming  that  differential  corrections  continue 
to  be  received  at  all  ranges  from  the  monitor.  Data  shown  in 
both  figures  are  for  a  composite  mechanization. 

Option  3  -  As  shown  in  Fig.  19.  application  of  informa 
tion  from  a  monitor  actually  degrades  performance  at  distances 
of  more  than  500  nm  from  the  monitor.  At  distances  in  excess 
of  500  nm,  best  performance  (shown  by  the  line  labeled  "no 
differential  corrections"  in  Fig.  19)  is  obtained  by  ignoring 
monitor  information  and  navigating  using  only  PPCs  for  the 
indicated  position. 

3.3.3  Summary 


Table  7  summarizes  the  recommended  techniques  for 
enroute  navigation.  A  relatively  high  degree  of  confidence 
can  be  placed  in  the  recommended  techniques  for  short  range, 
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RECOMMENDED  ENROL 
NAVI CAT I 


T3^^1 


DISTA.NCE  FROM 
•MONITOR 

DIFFERENTIAL 

CORRECTION 

RECEIVABLE"’ 

RECO.MMENDED 

OPERATIONAL 

TECHNIQUE 

<  1 2  5  nm 

Yes 

Apply  current  differential 
update 

125-500  r.m 

No 

.Modify  the  last  differ¬ 
ential  update  with  PPCs 
applicable  to  the  monitor 
site  and  the  current  air¬ 
craft  position  before 
application  to  the  meas¬ 
ured  phase 

>  500  nm 

No 

Ignore  the  monitor  infor¬ 
mation  and  navigate  using 
only  local  PPCs 

where  a  differential  correction  is  available,  and  long  range, 
where  monitor  information  is  not  applicable.  At  medium  range, 
however,  results  are  expected  to  be  sensitive  to  modeling 
errors,  especially  in  the  correlation  time  parameter  tj.  In 
addition,  it  has  been  assumed  that  PPCs  are  applied  with  no 
error.  It  is  recommended  ».hat  the  medium  range  navigation 
technique  be  tested  to  verify  the  predicted  performance.  (An 
algorithm  for  implementing  these  recommendations  is  described 
i n  Sect  ion  ^ . 3 . ) 
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IMPLEMENTATION  IN  AIJ\SKA 


Principal  requirements  for  implementing  a  differential 
Omega  navigation  system  in  Alaska  which  is  capable  of  supporting 
the  non-precision  approach  criteria  are  described  in  this  chap¬ 
ter.  Keconunendai ions  are  based  on  a  comparison  of  results 
presented  in  Chapter  3  with  specific  navigation  accuracy  re¬ 
quirements  described  in  Section  4.1. 


^4.1  ACCURACY  REQUIREMENTS 
4.1.1  Ove  rW  ew 

Navigation  accuracy  requirements  have  been  established 
by  the  FAA  for  aircraft  using  the  National  .Airspace  System 
(NAS).  The  accuracy  requirements  fulfill  two  purposes:  to 
ensure  safe  operation  of  individual  aircraft  and  to  provide 
control  of  the  airspace  to  afford  maximum-density,  safe  oper¬ 
ation  for  multiple  users. 

Two  distinct  categories  of  system  navigation  accuracy 
requirements  can  be  identified: 

•  Cert  1 f i cat  ion  requirements 

•  Airspace  planning  requirements. 


Each  category  has  its  own  particular  application. 

Certification  requirements  are  established  by  the  FAA 
in  order  to  qualify  equipment  for  particular  applications. 
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Certification  requirements  relate  to  equipment  and  represent 
the  most  stringent  category  ot  accuracy  requirements.  An 
example  of  certification  requirements  are  those  contained  in 
Ref.  8  for  use  in  certifying  Area  Navigation  Systems. 

.Airspace  planning  requirements  are  used  by  the  FAA  in 
establishing  procedures  for  aircraft  to  follow  in  particular 
circumstances;  e.g.,  in  preparing  a  Standard  Instrument  Depar¬ 
ture  \SID;  or  a  Standard  Terminal  Arrival  Route  (ST.AR).  These 
requirements  represent  planning  factors  associated  with  partic¬ 
ular  equipment  configurations  and  are  used  to  route  traffic 
around  obstacles,  define  .Minimum  Descent  Altitudes  (MDA),  and 
establish  navigation  fixes  (e.g.,  Initial  .Approach  Fix  (lAF) 
and  Final  .Approach  Fix  (FaF)).  Planning  requirements  are  used 
in  documents  like  Kef.  6  to  aid  in  designins  airspace  control 
procedures . 

Both  types  of  navigation  accuracy  requirements  are 
pertinent  to  the  discussion  of  the  Al aska/'Vukon  differential 
Omega  experiment.  If  a  differential  Omega  system  does  not 
meet  certification  requirements  under  the  Area  Navigation 
category,  the  system  cannot  be  used  as  a  primary  means  of 
navigation.  Regardless  of  what  accuracy  can  be  achieved,  it 
must  be  quantified  in  terms  of  planning  requirements  before 
procedures  can  be  developed  for  use  of  differential  Omega  in 
the  NAS. 

^.1.2  Non-Precision  Approach  Aid 

By  definition,  use  of  a  differential  Omega  navigation 
system  as  an  approach  aid  falls  in  the  non-precision  approach 
category;  no  provision  has  been  made  to  generate  an  electronic 
glide  slope.  Instrument  certification  requirements  for  non¬ 
precision  approach  from  Ref.  8  are  0.30  nm  along-track  and  0.3^ 


nm  cross-track  with  95%  probahility.  For  definiteness  in 
this  report,  these  accuracy  requirements  are  applied  for  an 
FAF  at  2  nm  from  the  airport  facility. 

Procedural  requirements  are  defined  by  system  capa¬ 
bilities.  Therefore,  the  data  presented  in  Chapter  3  of  this 
report  specify  the  procedural  accuracy  requirements  for  a 
differential  Omega  system.  For  comparison,  the  l.AF  for  VOR 
approach  to  Deadhorse  ,  .Alaska  corresponds  to  0.63  nm  (95%) 
along-track  and  0.50  nm  (95%)  cross-track  according  to  the  cri¬ 
teria  of  the  Kef.  6. 

4.1.3  Enroute  Aid 

The  instrument  certification  requirements  for  enroute 
navigation  from  Kef.  8  are  1 . 5  nm  (95%)  for  both  along-track 
and  cross- track  components.  Procedural  requirements  are  not 
applicable  in  this  category. 

4.1.^  Summary 

Two  types  of  navigation  accuracy  requirements  apply 
to  the  differential  Omega  e.vperiment:  certification  under 
.Area  Navigation  requirements  (Ref.  6)  and  procedural  require¬ 
ments  as  u.sed  in  the  planning  documents  (Kef.  6).  These  re¬ 
quirements  are  summarized  in  Table  8.  The  procedural  require¬ 
ments  data  1.S  an  example  obtained  by  analyzing  the  lAF  for  a 
VOR  approach  to  Deadhorse.  .Alaska. 


*The  allowance  for  Flight  Technical  Error  (PTE)  has  been 
removed  from  the  cross-track  error  component. 


TABLE  8 

navigation  accuracy  requirements  (95%  PROBABILITY) 


PROCEDURAL 

AID  C.ATEGORY 

COMPONE.NT 

CERTIFICATION 

EXAMPLE 

.Along-Track 
Error,  nm 

0.30 

0.63 

Non- Precision 

.Approach 

Cross-Track 
Error,  nm 

0.34 

0.50 

Along- t  rack 
Error,  nm 

1.50 

t 

NA 

Enrout  e 

Cross- 1  rack 
Error,  nm 

1.50 

t 

NA 

^NA  =  Not  Applicable. 


^.2  THE  ALASlOX  E.XPERIMENTAL  EQUIPME.N’T 

Three  types  of  equipment  are  needed  to  support  the 
use  of  differential  Omega  as  a  non - prec i s ion  approach  aid: 


•  Ground  station  equipment 

•  Data  link  equipment 

•  Airborne  navigation  equipment. 

Descriptions  of  each  of  these  equipment  categories  are  pre¬ 
sented  below. 

4.2.1  Ground  Station  Equipment 

Ground  station  equipment  must  include  Omega  monitor 
receivers,  data  processing  and  formatting  equipment.  The  ground 
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station  configuration  designtd  for  the  Al aska/Yukon  experiment 
IS  considered  adequate  to  support  certification  of  differential 
Omega  as  a  non-precision  approach  aid  if  it  is  increased  to 

ir 

provide  a  full  t i ghi - st at i on  capability  .  Nevertheless,  the 
analysis  presented  in  previous  chapters  validates  the  require¬ 
ment  for  certain  design  features  of  this  equipment  which  de¬ 
serve  special  emphasis: 


•  Precision  Icical  oscillator  at  the 
mi>nitor  (Rubidium  standard) 

•  Nearly  continuous  phase  and  amplitude 
data  (every  10  sec) 

•  Some  filtering  of  data  prior  to  trans¬ 
mission  (average  data  for  five  min). 

(.Actual  equipment  characteri  st  ics  are  shown  in  parentheses.) 

The  monitor  receivers  used  in  the  differential  Omega 
experiment  are  standard  TKACOR  599R  receivers  which  have  been 
modified  to  provide  a  digital  phase  output.  The  receivers  are 
augmented  by  a  precision,  external  atomic  frequency  standard. 
Therefore,  the  receivers  are  capable  of  measuring  single  sta¬ 
tion  Omega  phase,  and  independent  differential  phase  correc¬ 
tions  can  be  generated  for  each  Omega  signal. 

The  monitors  also  provide  nearly  continuous  phase  and 
amplitude  data.  The  data  are  filtered  to  smooth  out  noise  and 
then  sent  to  the  data  formatting  equipment.  The  effect  of 
this  filtering  is  depicted  in  Fig.  10. 

Signal  amplitude  data  are  used  in  the  ground  station 
equipment  to  derive  a  signal -qual i ty  index.  This  index  allows 
the  airborne  equipment  to  determine  the  quality  of  differential 
corrections  and  avoid  processing  signals  which  are  unreliable. 


♦The  experiment  will  use  only  signals  from  Norway,  Hawaii, 
North  Dakota  and  Japan. 


Differential  corrections  are  computed  by  subtracting 
the  received  Omega  phase  from  the  nominal  phase  for  the  moni¬ 
tor  site.  The  correction  is  then  formatted  and  inserted  an 
the  data  message.  .New  corrections  are  computed  every  10  sec. 

^.2.2  Data  Link 

The  data  link  designed  for  the  differential  Omega 
experiment  in  Alaska  appears  to  be  adequate  to  support  certif¬ 
ication  requirements  for  a  non-precision  approach  aid.  The 
message  format  should  be  expanded  to  allow  correction  proces¬ 
sing  for  all  eight  Omega  transmitting  stations,  but  sufficient 
channel  capacity  exists  for  this  modification.  (The  design 
channel  capacity  provides  an  information  rate  of  60  characters 
per  sec .  ) 

*^.2.3  .'Nirborne  Nav iJl^^i^n  Equ i pment 

Performance  predictions  contained  in  Chapter  3  indi¬ 
cate  that  the  TRACOR  7620  receiver  may  meet  certification  re¬ 
quirements  as  an  enroute  aid  to  navigation  under  the  criteria 
of  Ref.  8.  Uncertainty  in  the  value  of  the  correlation  coeffi¬ 
cient,  0.  which  cannot  be  resolved  without  processing  data 
from  the  e.xperiment,  provides  a  range  of  accuracy  numbers; 

★ 

e.g.,  cross-track  position  error  near  Inuvik  is  1.2  nm  (95%) 
for  p=0.9  and  2.3  nm  (95%)  for  p=0.5  compared  to  a  requirement 
of  1.5  nm.  The  true  correlation  coefficient  is  thought  (see 
Appendix  B)  to  be  much  nearer  1.0  than  0.5,  and  p=0.9  may  pro¬ 
vide  reasonable  predictions  of  the  TRACOR  7620  navigation  set 
performance . 

♦The  results  in  Chapter  3  are  expressed  as  rms  errors.  For  a 
joint  Gaussian  distribution  with  zero-mean,  identically  distributed, 
independent  components,  95%  probability  corresponds  to  a  circle 
of  radius  2.95  times  the  standard  deviation,  o,  of  each  component. 
The  95%  region  for  a  single  Gaussian  distribution  is  il.96o. 


The  predictions  in  Chapter  3  indicate  that  the  accu¬ 
racy  obtainable  with  the  TKACOh  7620  receiver  in  the  differen¬ 
tial  Oniega  configuration  does  not  meet  certification  require¬ 
ments  as  a  non- prec 1 s ion  approach  aid  under  the  criteria  of 
Ffel.  8.  For  example,  the  cross-track  error  ranges  between 
0.31  and  0.49  nm  (95%)  as  compared  to  a  requirement  of  0.30  nm 
(95%).  .Vevertheless ,  two  alternate  implementations  are  avail¬ 
able  which  offer  the  potential  to  meet  certification  require¬ 
ments:  operation  with  13.6  kHz  signals  only,  and  an  "optimal” 

two- f requency  implementation.  (Both  of  these  mechanizations 
are  discussed  in  Section  4.3.) 

Regardless  of  how  the  navigation  receiver  is  mecha¬ 
nized,  It  must  provide  additional  features  to  the  pilot: 


•  Waypoint  flight  path  computation 

•  Guidance  display. 

The  waypoint  features  of  the  TKACOR  7620  set  are  typical  of  a 
large  number  of  commerc i al  mi  1 1 1 ary  receivers  and  provide  for 
operator  insertion  of  up  to  10  waypoints.  The  operator  can 
select  displays  which  provide  di st ance- to-go  .  desired-track- 
angle,  est imat ed- t ime-of - arr i val  ,  etc.  (To  be  most  useful. 
F.AA-deve loped  SID  and  STAR  procedures  would  include  tabulated 
coordinates  for  procedural  waypoints.)  These  displays  offer 
guidance  information  to  the  pilot,  but  more  standard  control 
features  are  desirable:  e.g.,  outputs  to  drive  a  course  devia¬ 
tion  indicator  or  radio  magnetic  indicator. 


4.3  OTHER  EXISTING  EQUIPMENT 

The  equipment  constructed  for  the  differential  Omega 
experiment  in  the  Alaska/Yukon  region  is  believed  to  be  the 
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only  equipment  currently  .ivujluble  for  use  in  an  airborne  ap¬ 
plication.  All  of  the  equipiiit-nt  is  special-purpose,  and  of  a 
developmental  nature. 

During  the  course  of  this  contract,  TASC  contacted 
numerous  manufacturers  of  Omega  navigation  equipment.  .Many  of 
them  expressed  an  interest  in  the  potential  market,  but  none 
disclosed  plans  to  develop  or  market  receiver  equipment  in  the 
near  future.  The  primary  concerns  expressed  by  equipment  manu¬ 
facturers  were  an  uncertainty  in  the  size  of  the  potential 
market,  the  absence  of  plans  and  standards  for  ground  station 
and  data  link  equipments,  and  the  marginal  nature  of  Omega 
coverage  in  large  sections  of  the  continental  United  States. 

At  the  same  time,  several  manufacturers  were  eager  to 
discuss  their  product  line  of  Omega  and  Omega/\’LF  communica¬ 
tions  receivers.  These  receivers  are  primarily  designed  for 
the  Commercial  market  on  long  over-water  flights  where  VOK/D.ME 
signals  are  unavailable. 


a. 4.  SUGGESTED  l.MPKOVEMENTS 

The  ground  station  and  data  link  equipments  described 
in  Section  .  1  are  adequate  to  support  analytical  certifica¬ 
tion  requirements  for  differential  Omega  operation  in  the  non- 
precision  approach  category  if  they  are  modified  to  provide  a 
full  e i ght - St  at  ion  capability.  Nevertheless,  experience  gained 
as  a  result  of  the  experiment  may  suggest  modifications  to 
improve  reliability  or  convenience. 

Two  alternative  approaches  to  the  mechanizat i on  of 
the  airborne  navigation  equipment  appear  to  offer  promise: 
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•  Single- frequency  operation  in  the  dif¬ 
ferential  Omega  mode 

•  Correction  including  predicted  propaga¬ 
tion  effects. 

Both  techniques  are  described  below. 

4.4.1  Single- Frequency  Operation 

The  composite  Omega  procedure  used  by  TR,^COK  in  the 
7620  set  IS  designed  to  reduce  receiver  sensitivity  to  certain 
types  of  propagation  anomalies  (e.g.,  Sudden  Ionospheric  Dis¬ 
turbances)  and  to  improve  performance  with  a  simplified  propa¬ 
gation  prediction  model.  Within  the  differential  Omega  region, 
the  ground  station  equipment  provides  a  more  accurate  means  of 
correcting  for  these  effects,  and  the  performance  degradation 
associated  with  processing  additional  independent  receiver 
noise  is  not  offset  by  a  corresponding  gain  in  navigation  per¬ 
formance.  (The  composite  mechanization  subtracts  the  phase 
received  on  13.6  kHz  from  that  received  at  10.2  kHz.  While 
this  reduces  propagation  errors,  it  increases  the  effect  of 
measurement  noise  because  the  noise  contributions  to  each  fre¬ 
quency  are  independent.) 

For  e.xample.  compare  the  position-fix  quality  in  terms 
of  cross-track  position  errors  for  the  TRACOR  composite  (p=0.9) 
and  the  s ingl e - f requency  (13.6  kHz)  mechanizations  in  Table  4. 
The  rms  random  propagation  errors  are  nearly  equal:  0.16  and 
0.14  nm .  respectively.  The  rms  receiver  errors  are  substan¬ 
tially  different:  0.22  and  0.07  nm ,  respectively.  The  total 
rms  errors  are  0.27  and  0.16  nm,  respectively.  This  difference 
is  directly  attributable  to  the  increased  sensitivity  of  the 
composite  mechanization  to  receiver  noise  as  compared  to  the 
single- frequency  mechanization. 
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The  use  ot  s  i  n^l  e- f  t  i-quency  ddtu  xn  ihe  differential 
Ofuega  region  improves  navigation  performance  for  the  non-pre- 
cision  approach  application,  l)ui  it  does  not  eliminate  the 
need  for  a  multiple  frequency  receiver  in  enroute  applications 
(Multiple  frequencies  allow  the  receiver  to  resolve  much  wider 
lane  ambiguities,  which  might  be  required  after  a  temporary 
in-flight  electrical  power  transient.)  .Multiple  frequencies 
allow  better  operation  outside  the  differential  Omega  region 
and  provide  potential  insensitivity  to  propagation  anomalies. 

u  .  I4 . 2  Propagation  Model  Correction 

If  the  airborne  receiver  is  supplied  the  coordinates 
of  the  ground  station  monitor,  a  propagation  model  correction 
algorithm  can  be  implemented  which  offers  the  potential  for 
significantly  improving  the  transition  from  differential  Omega 
to  Omega  navigation  outside  the  differential  Omega  region. 

This  IS  similar  to  the  procedure  described  in  Ref.  9,  and  it 
is  supported  by  the  results  given  in  Figs.  17  and  18. 

The  mechanization  of  a  propagation  model  correction 
algorithm  might  have  the  following  form.  Denote  by  .  the 
differential  correction  broadcast  by  the  ground  station  for 
station  1  on  frequency  j.  Then 


t  is  the  measurement  instant 

C-  (x)  is  the  nominal  Omega  phase  at  x  for  station  i 
^  "  at  frequency  j 

s'"  i(x>()  fs  the  measured  Omega  phase  at  x  at 
instant  t 
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Denote  the  receiver-computed  propagation  correction  (skywave 

correction  for  the  TKACOR  7(>20)  by  ♦P  (x.t)  and  define  the 

^  •  J 

propagation  model  correction  as  follows: 


^♦P  .(t  )  =  ♦P  .(x.t)  -  .(x.t  ) 


(  K) 


It  the  differential  corrections  are  lost,  perhaps  because  the 
aircraft  has  departed  the  differential  Omega  region,  define  a 


new 


propagation  correction  (v.A)  (for  location  y  and  time 

rger  than  t^.  the  tim< 
unava 1 1 ab 1 e )  according  to 


V  larger  than  t^.  the  time  at  which  new  corrections  become 


-^.(y.x)  =  e-^.(y.M  - 


(15) 


where  3  is  a  weighting  factor. 

A  candidate  algorithm  for  computing  3.  which  is  con¬ 
sistent  with  the  data  in  Chapter  3.  is 


6  = 


(16) 


where  T  is  a  memory  constant.  The  choice  of  T  depends  on  a 
m  ■  ro 

number  of  factors  including  the  Omega  characteristics  of  the 
operating  region  and  the  flight  characteristics  of  the  air¬ 
craft.  (A  reasonable  value  for  T  is  expected  to  be  45  min.) 

m 

The  algorithm  described  by  Eqs.  14  through  16  provides 
one  method  for  implementing  the  recommendations  of  Section  3.3. 
It  has  the  advantage  that  the  transition  between  applying  a 
phase  monitor  correction  and  operating  only  with  computed  propa¬ 
gation  corrections  is  smooth,  and  the  pilot  will  not  experience 
rapid  jumps  in  the  indicated  Omega  position.  The  potential 
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di sadvantages  of  the  algorithm  are  that  it  has  not  been  tested 
and  it  may  be  sensitive  to  assumptions  about  the  nominal  param¬ 
eter  values  in  the  Omega  phase  error  model. 


CONCLUSIONS 


3  . 


This  report  presents  the  results  of  an  analytic  eval¬ 
uation  of  the  differential  Omega  experiment  in  the  Alaska/ 

Yukon  region.  Primary  emphasis  is  placed  on  the  particular 
equipment,  installation  and  procedures  peculiar  to  this  experi¬ 
ment.  The  analysis  also  includes  an  as  .essment  of  the  expected 
performance  of  an  "optimal"  two- f requency  receiver  using  the 
experimental  differential  Omega  correction  format.  The  latter 
results  provide  an  analytical  justification  of  the  suitability 
of  differential  Omega  for  non-precision  approach  applications. 
This  analysis  is  based  on  a  worldwide  Omega  error  model  which 
has  not  been  validated  for  the  Alaska/Yukon  region. 

Performance  results  presented  in  Chapter  3  for  nomi¬ 
nal  parameters  of  the  Omega  error  model  described  in  Chapter 
2,  indicate  that  the  experimental  differential  Omega  equipment: 


•  .'fay  not  meet  certification  requirements 
as  an  Area  Navigation  System  (Ref.  8) 

for  non-precision  approach  without  further 
modi f 1  cat  ions 

•  Should  meet  certification  requirements 
as  an  Area  Navigator  for  enroute  navi¬ 
gations 

•  Will  not  be  adversely  affected  by  the 
loss  of  signals  from  any  one  of  four 
Onega  transmitting  stations  normally 
received  in  the  experimental  region 

•  Is  expected  to  maintain  adequate  Omega 
signal  tracking  during  standard  rate 
turns . 


APPENDIX  A 

OMEGA  PHASE  ERROR  MODEL 


A.l  OPTIM.AL  OMEGA  SIGNAL  PROCESSING 

A  model  is  developed  in  this  section  for  the  position 
error  resulting  from  optimal  (minimum-variance)  usage  of  phase 
information  from  two  Omega  frequencies,  10.2  kHz  and  13.6  kHz, 
for  each  of  four  stations; 


• 

A 

( Norway ) 

a 

C 

(Hawai i ) 

• 

D 

(North  Dakota) 

a 

H 

(Japan ) . 

Based  on  the  estimated  position  of  the  navigation  set.  an  es¬ 
timated  phase  is  computed  for  each  station  i  on  each  fre¬ 
quency  j  as  follows; 


(A-1) 


where 


R  IS  the  range  to  station  i  from  the  estimated 
^  position 

X.  is  the  nominal  wavelength  of  a  signal  with 
frequency  j 


The  difference  between  estimated  phase  and  measured  phase 

•  is  the  phase  deviation: 


Phase  deviation  is  related  to  position  error  as  follows 


■^^1  .10.2 

COSOj 

^10.2 

s  1  n4i  ^ 

^10.2 

-  “ 

*^1,10.2 

•^®2.10.2 

cos(;i2 

siniii., 

*^2,10.2 

^10.2 

^10.2 

•^*‘’3.10.2 

COSlP^ 

^10.2 

Si  nui^ 

^10.2 

*®3.10.2 

•^^A.10.2 

i  -^*1.13.6 

! 

cosw^ 

^10.2 

COSVP  J 

^13.6 

s  i  niii^ 

^10.2 

s  i  nv  2 
^13.6 

AE 

^'^A.10.2 

*^1,13.6 

j 

1  -^^2.13.6 

1 

cost*' 

SI  nu>2 

*^2,13.6 

^13.6 

^13.6 

•^®3.13.6 

cosiii  2 

^13.6 

sintp-j 

^13.  t> 

**3,13.6 

•^*4.13.6 

cosv^ 

^13.6 

s  1  nv^ 

^13.6 

m 

**A.13.6 

where 


•^1 

AN 

AE 


is  the  phase  deviation  from  station  i  on 
frequency  j 

is  the  bearing  to  transmitter  i 
IS  north  position  erior 
is  east  position  error 


ij 


(A- 


is  the  phase  measurement  error  from  station 
1  on  frequency  j. 


Phase  measurement  error  is  comprised  of  phase  error  from 

all  sources  except  position  error,  and  includes  receiver  error, 
clock  synchronization  and  drift  error,  and  signal  propagation 
prediction  error.  .Note  that  Eq .  A~3  applies  in  the  non-differ* 
ential  mode  as  well  as  the  differential  mode.  In  the  differen¬ 
tial  mode,  a  correction  from  a  nearby  monitor  station  is  applied 
to  the  measured  phase  0*1*^  in  .m  attempt  to  minimize  the  phase 
measurement  er-or  It'  *  he  <’i  f  ferent  i  al  mode,  repre¬ 

sents  the  residual  phase  measuretent  error  after  correction  by 
the  monitor 


Equation  .^-3  can  be  rewritten  as  follows; 


=  H  Ax  ♦ 


It  can  be  shown  (het.  10)  that  the  minimum-variance  estimate 
of  given  the  phase  deviations  ^  in  the  presence  of  phase 

measurement  errors  is 

^  =  (H^R‘^H)'^H^R*^Ad  (.A-5) 

\  T 

where  R  -  E(^  6o  )  and  £(•}  is  the  expectation  operator.  The 
covariance  of  the  estimate  ^  is  computed  as  follows: 

P  =  (H^R'^H)'^  (A-6) 

where 

P  El  (^  -  ^)(^  - 

To  evaluate  the  variance  in  the  estimate  of  position  error  an 
expression  for  the  matrix  R  is  required. 


A-3 


To  determine  R,  consider  the  phase  measurement 

error  from  station  i  on  frequency  j  caused  by  all  sources  ex¬ 
cept  position  error.  .Assume  that  the  navigator  is  at  position 
and  obtains  a  phase  measurement  at  time  t  .  Then  the  meas¬ 
ured  phase  can  be  expressed  by 


«^‘j(Xj  .tj  )  = 


or 


(A-7) 


1 J 


ur.e  re 


(A-8) 


IT 

C  ‘  is  the  measured  phase  from  station  i  on 
frequency  j 


0 


1  j 


is  the  "nominal”  phase  determined  from  range 
to  transmitter  and  signal  wavelength 


6«  is  phase  measurement  error  caused  by  all 
sources  except  position  error 

6c^  is  the  phase  perturbation  caused  by  predictable 
propagation  effects 


6o 


r 

ij 


is  the  phase  error  caused  by  unpredictable 
(random)  propagation  error 


is  the  constant  offset  portion  of  measured 
phase  introduced  by  the  receiver  (i.e., 
clock  error) 


6®'.  is  the  variable  offset  portion  of  the 

^  measured  phase  introduced  by  the  receiver 
(receiver  noise) 


Assume  that  a  monitor  measures  phase  at  position  X2  at  time 
t.).  The  measured  phase  at  the  monitor  is  given  by 


(A-9) 


*  6^,;(x-j,l-. )  ■*  60-.(x<jit-j) 

ij  -*.  ij  _ 

The  monitor  correction,  transmitted  to  the  navigator,  is  roeas' 
ured  phase  minus  nominal  phase: 


CORRECTION  =  ^"'.(x,,t^)  -  ®,.(x.) 

1  j  —z  k  1  j  ^ 

=  6«>P.(Xp,tT)  •»  6«i^j(X2.t2)  6ci^.(x^,t2) 


^0^j(x2 ,t2) 


(A-10) 


The  navigation  set  receives  the  correction  and  uses  it  to  esti 
mate  nominal  phase  as  follows: 


«ij(xj.ti)  =  ^"y-1’^1^  * 


(A-11  ) 


Substituting  Eq .  A-10  into  Eq .  A-ll  yields 


ij'i2 


(A-12) 


Combining  Eqs .  A-12  and  A-8  yields  an  expression  for  the  error 
in  measured  phase  at  the  receiver: 


^♦^j(Xi  .tj) 


♦ij(Xi.ti)  -  ♦ij(Xj) 


=  *  6»Pj(X2,t2) 
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*  *^1  ^ 


-  -  6»^.(X2.l2) 


(A-13) 


Equation  A-13  shows  that  measured  phase  error  depends 
only  on  the  di f  ference  in  the  four  types  of  errors  at  the  navi¬ 
gator  and  the  monitor.  The  monitor  is  effective  because  the 
terms  *»tid  change  relatively  little  over  short  dis¬ 

tances  and  short  time  intervals;  i.e.,  these  errors  are  highly 
correlated  in  time  and  space.  The  other  terms,  and 

which  are  uncorrelated  between  the  monitor  and  the  navigator, 
must  be  relatively  small  for  the  monitor  to  be  effective. 
Equation  A-13  is  used  to  compute  the  matrix  R. 


The  following  assumpt i ons  simplify  the  analysis: 

1.  fto^j(XTt2)  =  PPCs  are 

assumed  to  not  change  over  the  distance 
j,\2*Xi  1  and  the  time  1^2’^ir 

2.  6c^j(x2.t,)  -  =  0-  Any  bias 

offset  in  phase  introduced  by  the  re¬ 
ceiver  IS  assumed  to  be  negligible. 

3.  Random  propagation 
anomalies  are  assumed  to  be  uncorrelated 
with  receiver  noise. 

U,  E16c  .  =  0.  Any  two  errors  involving 

1  J  K  J 

two  different  transmitters  are  assumed  to 
be  uncorrelated. 


Under  these  assumptions,  it  can  be  shown  that  the  matrix  R  is 

2 

comprised  of  common- frequency  terms,  of  the  form  E{|6*^jJ  J, 
and  cross- frequency  terms,  of  the  form 


A-6 


Conuaon-  f  rfequency  terms  are  in  turn  composed  of  two 
types  of  terms:  receiver  error  and  random  propagation  error. 
The  model  for  variable  receiver  error  is  chosen  to  be  a  dis¬ 
crete  while  noise  sequence  with  a  variance  of 

El  I6«^^'.(x.t  )  1^1  =  aj  (A-1^) 

The  assumed  model  for  random  propagation  error  is  adapted  from 
Ref.  6  and  applies  to  the  case  ^2*^1  2  hr  for  a  frequency 

of  10.2  kHz: 

I*®  i  .  10. 2^-2  •  ^  .  10. 2^-1  ^  ^  ^ 


(A-15) 

where 


d  = 

1500  nm 

9.83 

cec 

^  = 

0.943  hr 

17.8 

cec 

T  — 

•  1  • 

12  hr 

8.13 

cec 

= 

8  hr 

1.98 

cec 

The  model  is  calibrated  with  data  at  10.2  kHz.  At  13.6  kHz, 
It  is  assumed  that  the  same  functional  form  applies  but  in¬ 
cludes  a  weighting  factor.  V: 
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■J  /  liv'-lM^  2  /I* 

=  W‘  exp  (-  ■  -■d  ■  h^i 


^  0,^  -  03-  cos  It2  -  ij 


„(!i-.J) 

■1')"  ‘’4^  <^‘>5  ( 


COS  ^  U 

\  '4 


(A-lb) 


2  -  '1')' 


Cross- f requency  terras  involve  only  random  propagation 
error.  From  Eqs .  A- 15  and  A-lb.  it  can  be  shown  that  the  cross- 
trequency  term  js  given  by: 


Ell  (be^_j(x2.t2)l  l6o[j^(Xj  .tj))) 

I  I  ^2  '  -1 1  2  /  *  ^2  '  1  ^ 

=  pv  exp  (- 

”)  2  /  ’’n  \ 

"  £^2^  "  °3  '^2  *  '1') 


>  0^2  |2£ 


for  j?'k 


(A-17) 


where  p  is  the  correlation  coefficient  between  frequencies. 


The  matrix  K  is  composed  of  common-  and  cross-fre¬ 
quency  terms  as  follows: 


0 


0 


0 


0 


0 


where 


kj  =  El  |6*"Jj(X2.i2))^}  ♦  El  .t  j  )  r  } 

♦  2EI 

-  2E{ I A« ^  jQ  2^52  •  ^2 ^ ,10.2^-1 ’^1 ^ 

Ro  =  El l6oyj(x2.t2)l’)  ♦  Ej |fi*jj(Xj  .tj  )  1^) 

•  2EII*<.:3.6<£.>'I^I 

2E|  |6*  J  13  ^(X2  .  t2  )  I  l**i  ,i3.(,^5i  **^1  H  1 
R3  =  2E|  |6«[j(x.t  )  1  I6*[j^(x.t)  1) 

-  2E|  16c^(x2.t2))  I6^^^(xj  .tj  )  J1  (ji^k) 

A. 2  TRACOR  COMPOSITE  SIGNAL  PROCESSING 

In  the  previous  section,  the  covariance  of  the  posi¬ 
tion  estimation  error  was  derived  (Eq.  A-6)  assuming  optimal 
use  of  eight  phase  measurements  (Eq.  A-3).  In  this  section, 
the  covariance  of  the  position  estimation  error  is  calculated 
assuming  the  eight  phase  measurements  to  be  combined  as  in  the 
TRACOR  7620  receiver  to  form  three  "composite”  phases  of  the 
form 


^1 j  ,c  ”  ^^*i  ,13.6 


■  ♦j,13.6^  ■  ^^♦i,10.2  ■  *3,10.2^ 


(A-18) 


where 
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i  j  .c 


IS  I  he  coniposiie  phase  from  I  ransmi  1 1  ers 
i  and  j 


1  . 1 3  .  f 


IS  the  phase  from  transmitter  i  on  fre 
quency  13.6  kHz 


.10.2 


IS  the  phase  from  transmitter  j  on  fre¬ 
quency  10.2  kHz 


a.b  are  constants  (a  =  15/8,  b  =  3/2) 


In  analogy  to  Section  A.l,  composite  phase  deviation, 
IS  related  to  position  error  as  follows 


ij  .c 


12. c 


.io, 


3  .c 


.io 
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cost*!  -C0SV2 


cose  2 -COSVj 


COSw  2  -  COSv»^ 
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12. c 


*♦13. c 


14  ,c 
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where 


^*ij.c  ^  '‘^*®i.l3.6  ■  *  ^^**i.l0.2  ’  **j.l0.2 


13.6  ^10 


1) 


-1 


=  composite  wavelength 


Phase  measurement  error  6*^^  (j  =  13.6  or  10.2)  is  given  in 
Section  A.l  (Eq.  A-13). 

Equation  A-19  can  be  written  in  the  form 


^  =  He  ^  ^c 


(A-20) 


A-IO 


The  minimum  variance  estimate  oi  Ax  is  then  given  by 


where 

«c  =  ^l*lc 

The  covariance  of  which  represents  a  lower  bound  for  the 
covariance  of  position  error  resulting  from  a  single  set  of 
three  composite  phase  measurements  using  the  TRACOR  7620  algo¬ 
rithm,  IS  given  by: 

Pc  =  (A-22) 

The  matrix  R  can  be  determined  in  a  manner  similar  to  the 
c 

method  in  Section  A.l  and  has  two  types  of  terms: 

•  Diagonal  .  E| (6# .  .  )^} 

*  J  • 

•  Off-diagonal.  E{16».j^^)  ^ 

Diagonal  terms  are  given  by: 

E{16»i.  cl^l  =  E{|6*[  jQ 

-  ^b  E{16^^  jQ  1^^^  jQ  2^—1'^!^^^ 

■  ,13.6*-2  •'^2*'  ,13.6*Sl ''l  ^ 

“  8ab  E|(6#j  22  10,2^  —  *^^^^ 
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-  Sab  El  |6^[j(x2.i2)) 

^  2(3^  b-)  E{  16»^’.(X2.t2)l^) 


^  2(a2  .  b^  El  |6*y  (xj.tj))^) 


• •  ' ij  -1 

for  (j  ^  k) 


(A-23) 


Terms  on  the  right  side  of  Eq .  A-23  are  given  by  Eqs.  A-14, 
A-15,  A-16,  a-17.  Of f-di agonal  terms  are  equal  to  one-half 
the  diagonal  terms: 


EII*«ij  c*  •^•ik.c^*  ■  2  '**ij  •  J  ^  ^  (A-24) 


The  matrix  is  written  as  follows; 


=  R 


where 
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APPENDIX  B 

FREQUENCY  CORRELATION  OF  OMEGA  ERRORS 


This  appendix  presents  the  approach  used  to  provide  a 
statistical  prediction  of  the  c ross- f requency  correlation  between 
the  propagation  medi um* i nduced  random  phases  of  10.2  and  13.6 
Omega  signals.  Both  the  10.2  and  13.6  kHz  signals  are 
assumed  to  be  radiated  by  the  same  Omega  transmitter  and  received 
simultaneously  at  the  same  measurement  site.  The  equation 
for  the  c ross - f requency  correlation  coefficient  is 


‘^'10.2.13.6  ■ 


_ ^|ft»10  2  ^^13.6* 

E(  (6*10.2^^' 


(B-1  ) 


where 


T  IS  the  propagation-medium*  induced  random 
phase  in  the  10.2  kHz  signals 

propgat  i on -medi um- i nduced  random 
phase  in  the  13.6  kHz  signals. 

At  a  given  receiver  location,  the  phase  of  a  received 
Omega  (or  VLF)  signal  is  a  function  of  signal  frequency,  signal 
(i.e.,  t ransmi 1 1 er- to- rece i ver )  path  length,  and  propagation 
medium  properties  along  this  path  such  as  ground  conductivity, 
geomagnetic  latitude,  bearing  angle  and  the  (D-region)  iono¬ 
sphere.  The  ionosphere  properties  vary  significantly  with  the 
solar  zenith  angle  along  the  path.  The  primary  source  of  propa¬ 
gation-medium-  induced  random  signal  phase  is  the  random  temporal 


♦The  one-second  delay  between  successive  Omega  bursts  is 
negligible  in  this  context. 
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variations  in  the  properties  of  the  ionosphere  along  the  path. 
Therefore,  the  subsequent  random  phase  analysis  considers  only 
random  changes  in  the  ionosphere  properties. 

At  VLF ,  the  D- region  ionosphere  properties  are  ade¬ 
quately  character ized ,  according  to  Wait  and  Spies  (Kef.  11), 
by  two  parameters:  the  ionosphere  reflection  layer  height,  h, 
and  the  ionospheric  conductivity  gradient,  p.  According  to 
Kelly  (Kef.  12),  each  parameter  can  be  assumed  to  be  independ¬ 
ently  variable  and  characterized  by  a  Gaussian  probability 
distribution.  “The  frequently-used  values  for  the  mean  and 
standard  deviation  of  the  ionospheric  parameters  are  given  in 
Table  B-1  (Kefs.  12  and  13).  for  both  day  paths  (i.e,,  local- 
noon  conditions  along  the  entire  path)  and  night  paths  (com¬ 
plete  darkness  along  the  entire  path).  The  mean  values  cor¬ 
respond  to  .1  "quiet"  ionosphere. 

Sample  phases  ftCj-  (f  =  10.2  or  13.6  kHz),  required 
for  computing  the  frequency  cross-correlation,  Eq .  B-l,  are 
obtained  by  employing  the  commonly- used  deterministic  VLF  sig¬ 
nal  propagation  prediction  model,  known  as  IPP  (Refs.  14  through 
16).  In  these  computations,  the  non-random  part  of  the  phase 


TABLE  B-l 

STATISTICS  FOR  IONOSPHERIC  MODEL  PARAMETERS 


PARAMETER 

HEIGHT 
h  (km) 

CONDUCTIVITY 

GRADIENT 

P(kro)'^ 

SOLAR 

ILLUMINATION 

MEAN 

f, 

STANDARD 

deviation 

°h 

MEAN 

B 

STANDARD 

DEVIATION 

Day 

72 

2 

0.3 

0.05 

Night 

87 

2 

0.5 

0.1 

is  associated  with  ’'iioniuial”  propagation  properties  ol  the 
signal  path  assuming  a  quiet  ionosphere.  The  random  part  of 
the  phase,  on  the  other  hatid .  is  assumed  to  be  induced  by  spa¬ 
tially  uniform,  random  changes  in  the  quiet  ionosphere. 

The  linearized  e.vpression  for  the  random  phase  due  to 
ar.  assumed  random  change  of  one  standard  deviation  in  the  iono 
spheric  parameters  h  and  p  is 

Be  f  9®  £ 

=  TF  °h  ^  “9P 

wnt  re 

h  =  f.  . 

6  =  P  •  Ojj 

The  phase  of  a  propagating  signal  is  approximated  by 


(B-2) 


(B-3) 


B-3 


(B-5) 


|af(Oh.Op 

)l  d 

where 

aVf 

^  '’h  ^ 

avf 

ap 

a ( argA  ^  ) 

/I  ^ 

3 ( argA  ^ ) 

ah 

‘^h 

ap 

No:^  Lq .  B-5 

is  derived  by 

assuming  constant 

(B*e>) 


(B-7) 


:orn  rrtridom  ionospheric  parameter  changes  <or  irregularities) 
.ilong  the  entire  signal  path  length  d.  If.  however,  the  iono¬ 
spheric  irregularity  does  not  extend  over  the  entire  signal 
path,  d  in  Eq .  B-^  should  be  replaced  by  the  actual  irregular¬ 
ity  length.  Furthermore,  the  excitation  factor  phase  contribu¬ 
tion  .  bj ,  should  be  (1)  multiplied  by  0.5  if  the  irregularity 
covers  either  the  transmitter  or  receiver  alone,  or  (2)  zero 
if  the  irregularity  covers  neither  the  transmitter  nor  the 
receiver.  Assuming  that  the  ionospheric  irregularity  covers 
the  entire  path,  the  equation  for  the  cross- frequency  correla¬ 
tion  coefficient  becomes 


10.2,13.6  = 


^^10.2  ^13.6^  *  ^^^10.2  ^13.6^ 

— T - T - -5 - T77 — T - T - >5 - - 

^^10.2  *  ^10.2^  "^^13.6  *  ^13.6^ 


1/2 


(B-8) 


The  correlation  coefficient  has  been  computed  for 
both  day  and  night  conditions  over  the  four  different  propaga¬ 
tion  paths  identified  in  Table  B-2.  These  paths  were  speci¬ 
fically  chosen  to  characterize  the  range  of  propagation  path 
properties  expected  along  the  Omega  signal  paths  of  interest 
in  the  Alaska/Yukon  region. 
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I  ABLE  B-2 

HKOEAf.AllON  CATH  PKOFEK  I  I  ES 


Table  B-3  presents  the  random  phase  terms,  a^  (Eq. 

B-bJ  and  b^-  ( Eq .  B-7)  corresponding  to  assumed  random  changes 
(lo)  in  the  ionospheric  parameters  (given  in  Table  B-1)  for 
each  of  the  tour  propagation  paths  identified  in  Table  B-2  at 
10.2  and  13. b  kli?. .  The  results  in  Table  B-3  indicate  that  the 
sign  of  either  coefticienl,  a^  or  b^ ,  is  not  altered  when  f 
changes  from  10.2  to  13.6  kHz.  This  means  that  both  terms  in 
the  numerator  of  the  correlation  coefficient.  Eq .  B-8,  are 
pc'Sitive.  and  the  correlation  coefficient  is  positive.  Further 
more,  since  b^  is  much  less  than  a^d  for  d  greater  than  100 
km.  the  results  show  that  for  signal  length  d  greater  than  100 
Km.  the  random  phase  can  be  approximated  by 


~  a^  d  (B-9) 


TABLE  B-3 

COMPETED  VALUES  OF  COEFFICIENTS  AND  bj- 
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Thus,  for  signal  path  lengths  greater  than  100  kni ,  the  cor¬ 
relation  coefficient  is  approximately  equal  to 

^10.2.13.6  :  tlL-lJjJ =  I  (B-10) 

^10.3  ^13.6 

Both  Eqs.  B*9  and  B-10  also  hold  even  if  the  irregularity  is 
smaller  than  the  signal  path  length,  but  both  the  signal  path 
and  irregularity  are  longer  than  100  km. 

In  conclusion,  the  analysis  presented  herein  indi¬ 
cates  that  the  c ross - f requency  correlation  between  10.2  and 
13. B  kHz  signal  random  phases  in  Alaska  is  positive  and  ap¬ 
proaches  unity  if  both  the  signal  path  length  and  the  random 
ionosphere  irregularity  length  are  greater  than  100  km. 
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APPENDIX  C 

TRACOR  7620  NAVIGATION  SET  PERFORM.ANCE  ANALYSIS 


C.:  TR.ACOR  FILTER 

C . 1  1  In: roduct ion 

The  TR.ACOK  Omega  navigation  algorithm  (Kef.  17)  em¬ 
ploys  a  di sc  re t e - 1 ime  position  filter.  Using  a  true  airspeed 
indicator  and  true  heading  information,  the  algorithm  propa¬ 
gates  a  dead- reckoned  position  on  a  1.25  sec  cycle.  On  a  10-sec 
cycle.  Omega  signals  are  used  to  provide  external  updates  to 
the  system.  The  mechanization  is  similar  to  that  of  a  Kalman 
filter  in  that  it  is  a  recursive  estimator.  However,  it  is 
suboptimal,  since  the  update  gam  matrix  is  not  dependent  on 
the  computed  current  error  covariance. 

Error  equations  are  presented  for  the  TRACOR  7620 
navigation  filter.  By  transformation  into  the  2  domain  (Ref. 
19).  certain  characteristics  of  the  system  may  be  examined 
analytical ly .  Among  these  are  system  stability,  damping,  and 
state  variable  coupling.  This  transform  approach  also  permits 
analysis  of  the  statistical  response  of  the  system  to  Omega 
phase  measurement  errors.  The  power  spectral  density  matrix, 
which  is  the  transform  of  the  error  correlation  matrix,  is 
presented  analyt i cal ly . 

C.1.2  Filter  Error  Equations 

The  filter  equations  may  be  written  in  10-sec  differ¬ 
ence  form  as 


x(kT 


1 

0 

10 

0 

cos(TH(kT-T)) 

0 

1 

0 

10 

x(kT-T)  - 

e 

sin(TmjkT-Tj))^ 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

TAS(kT-T) 


(C-1 ) 


and 

x(kT*)  =  x(kT’)  K(z  •  l^-(x(<T’))l  (C-2) 


whtre  kT  and  kT^  are  the  time  immediately  before  and  after 
updates  are  incorporated. 


The  state  estimate  is 

X  =  ILAT.  LON.  WlNDj^. .  WINDj.|^  (C-3) 

in  radians  and  rad i ans.  sec ,  as  appropriate,  and 

2^  (x)  is  the  vector  of  nominal  Omega 
composite  predictions  in  cycles 

z  is  the  vector  of  Omega 
composite  measurements 

V  is  the  vector  of  Omega  phase  measure- 
~  ment  errors 

TH  is  true  heading 

TAS  is  true  airspeed  in  knots 

R  IS  the  radius  of  the  earth  in  nautical 
^  miles 


T  is  the  time  inter\’al  ,  10  sec. 


The  form  of  the  gain  matrix  K  is 

K  =  Pf  H(x)R^*^  (C-C) 

where  and  Rf'^  art  fixed  matrices.  The  measurement  matrix 
H  can  be  used  to  approximate  the  difference  between  Omega  pre¬ 
dictions  and  measurements: 

2  -  l^Cx)  =  H(x)  lx  -  x)  (C-5) 

including  Andover -Lambert  ellipticity  corrections.  (For  sim¬ 
plicity,  the  functional  notation  H(x)  is  dropped  hereafter  and 
only  H  is  used .  ) 

By  employing  direct  perturbation  and  combining  the 
above  equations,  it  can  be  shown  that  the  error  equations  are 

^(kT)  =  II  -  KH)  l^^(kT-T)  Llw(kT-T)  -w(kT)l 

*  G  V  (kT))  ♦  K  v(kT)  (C-6) 

— m  — 

where  the  true  wind  vector  is 

w  =  iVINDj.,  VINDgl^  (C-7) 

the  instrument  error  vector  is 

V  =  UTAS.  .iTH)^  (C-8) 

'fli 

Omega  measurement  error  is  described  by  v  and  the  transition 
matrix  is 


where 


(C-9) 


f  =  -10  TaS  sin(TH)  scc(LAT)  tan ( LAT )/3600R^  (C-10) 


10  0  1  0 
0  10  0  1 


(C-11  ) 


cos^TH) 

sin(TH)  sec(LAT) 


-TAS  Sin(TH) 

TAS  cos(TH)  sec(LAT) 


(C-12) 

Equations  C-1  through  C-5  are  used  in  a  direct  simula¬ 
tion  of  the  TRACOR  7620  navigation  algorithm.  Equation  C-6 
provides  the  error  formulation  required  to  perform  statistical 
analysis  of  the  filter. 

C.1.3  Transform  Analysis 

If  the  measurement  matrix  H  is  held  constant,  Eq.  C-6 
can  be  Z-transformed  (Ref.  19)  and  manipulated  to  yield 

AX(i)  '  U1  -  (1  •  KH)*1‘^  1(1  -  z)(l  -  KH)  L  W(z) 


♦  z(l  -  KH)  G  V'  (z)  ♦  z  KV’(z) 
—re  ~ 


(C-13) 
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Although  H  is  geome try- dependent ,  it  tends  to  change  slowly  at 
moderate  aircraft  speeds.  Consequently,  this  transform  of  the 
"snapshot”  gives  a  good  representation  of  error  behavior,  and 
should  yield  considerable  insight  into  system  behavior  for 
other  flight  profiles. 

The  system  characteristic  function  was  calculated  for 
Inuvik,  assuming  availability  of  signals  from  stations  A,  C,  D 
and  H;  the  result  is 

c(z)  =  -  3.17z^  ♦  (3.6^  0.0^f)z^  -  (1.76  ♦  0.04f)z 


♦  (0.29  ♦  0.04f) 


(C-K) 


At  these  latitudes 

I f  |<  0.001 

For  f  =  0,  system  poles  are  found  at 
Zj  -  0.35 
Z2  =  0.96 

Z3  ^  =  0.93  1  j  0.027 

From  this  analysis  it  is  possible  to  conclude 

that 

•  The  system  is  stable 

•  The  ^2*^  pole-pair  gives  rise  to  a 
lightly  damped  oscillation 

•  Maneuvers,  which  change  the  value  of  f, 
will  not  affect  system  stability 

•  Manuevers  will  move  the  poles,  creating 
transient  errors. 
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C.1.4  Power  Spectral  Pgnsity  Analysis 


Given  a  transfer  function  and  power  spectral  densi¬ 
ties  of  input  disturbances,  it  is  possible  to  calculate  the 
power  spectral  density  (PSD)  matrix  of  the  error  vector  if  the 
processes  are  stationary  (Ref.  20). 

Under  this  assumption,  the  PSD  matrix  of  state  error  is 


P  (s) 

XX 


where 


IzI  -  (1-KH)  (  Iz-l 

-  Izl^  (1-KH)  GP  (s)5(  1-KH) 
m 


(1-KH)  LP^.(s)  L  (1-KH) 
IzJ^KR(s)R)  * 


|z  1  -  (1  -  KH)*J‘^ 


(C-15) 


z 

s 

Pv,.(s) 

tu 

R(s) 

(■) 


is  the  Laplace  transform  variable 

is  the  wind  PSD  matrix 

is  the  instrument  error  PSD  matrix 

is  the  Omega  measurement  error 
PSD  matrix 

indicates  the  complex  conjugate  of 
a  transposed  matrix 


Pj^^  is  a  steady-state  power  spectral  density  and  does  not 
reflect  the  lime-varying  error  covariance  propagation  charac¬ 
teristics. 


C.2  TRACOR  FILTER  WIND  RESPONSE 
C.2.1  Introduction 

One  factor  contributing  to  position  error  is  wind. 
The  impact  of  wind  on  position  error  covariance  is  difficult 
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to  decernine.  since  relatively  little  is  known  about  wind  cor¬ 
relations  in  the  Alaska/Yukon  region. 

Since  the  TRACOR  filter  carries  state  variables  for 
north  and  east  wind,  a  constant  wind  would  be  accurately  esti¬ 
mated.  but  changing  winds  will  induce  transient  errors.  Since 
the  only  observations  of  wind  are  from  induced  Omega  position 
error,  these  transients  may  decay  slowly. 

C.2.2  Winds  in  Alaska 

Data  were  analyzed  for  surface  winds  at  Eielson  AFB . 
near  Fairbanks.  Alaska  (Ref.  IS);  the  following  were  noted: 

•  The  air  was  calm  for  50%  of  the  observ'a- 
tions  on  a  year-round  basis 

•  Wind  speed  was  less  than  10  kt  for  97% 
of  observations 

•  No  observation  was  made  of  windspeed 
greater  than  3^  kt 

•  Mean  scalar  wind  speed  was  2.4  kt . 

However,  the  relation  of  surface  winds  to  wind  in  the 
upper  air  is  not  available,  and  the  time  correlation  of  wind 
speeds  in  different  directions  is  unknown.  As  a  result,  the 
contribution  of  wind  to  position  error  covariance  is  difficult 
to  determine. 

C . 2 . 3  Simulation  of  Wind  Effects 

In  order  to  assess  the  impact  of  wind  on  position 
error,  a  number  of  simulations  were  run.  The  simulations 
whose  results  results  are  depicted  in  Figs.  C-1  and  C-2  used 
wind  speeds  of  10  kt  and  50  kt,  respectively.  Other  simula¬ 
tion  conditions  were: 


C-7 


Initial  position:  68®N ,  134®W 
Zero  initial  wind 
True  airspeed:  250  kt 


•  Signals  from  transmitting  stations  A,  C, 

D,  and  K  assumed  to  be  available 

•  Perfect  reception  and  monitor  correc¬ 
tions  assumed 

•  Initial  track:  180® 

•  At  2  min,  north  wind  undergoes  a  step 
change  from  zero  to  the  appropriate  mag¬ 
nitude 

•  At  20  min,  wind  direction  shifts  to  east 

•  Also  at  20  min,  the  aircraft  starts  to 
turn  clockwise  through  270®  at  3®/sec 

Figure  C-1  and  C-2  show  the  same  features  predicted 
by  the  transform  analysis.  The  lightly  damped  oscillation  is 
not  a  dominant  feature  of  the  system  response.  Some  coupling 
exists  between  north  and  east  position  errors.  Also,  note 
that  maximum  radial  position  errors  occur  at  A : 50  and  22:A0 
flight  time  in  both  simulations  -  2:50  and  2:40  after  wind 
shifts . 

The  difference  between  the  results  of  the  two  simula¬ 
tions  is  in  error  magnitude.  The  maximum  position  errors  in¬ 
duced  by  the  wind  errors  are  provided  in  Table  C-1.  Note  that 
for  this  range  of  wind  magnitudes,  the  TRACOR  7620  response  to 
wind  shifts  is  nearly  linear. 

Position  errors  induced  by  wind  shift  are  important 
but  not  critical.  Maximum  position  errors  of  approximately 
0.23  nm  per  10  kt  of  wind  step  may  be  expected.  Should  wind 
change  mote  gradually,  smaller  maximum  errors  would  be  expec- 


*(min:sec)  of  simulation  time. 
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Figure  C-1  TRACOR  7620  Response  to  10  kt 
Wind  Shift 
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Figure  C-2  TRACOR  7620  Response  to  50  kt 
Wind  Shift 


ted.  While  slower  than  for  maneuver- induced  error,  transient 
decay  still  reduces  position  error  to  10%  of  maximum  within  10 
min  of  the  wind  shift. 


TABLE  C-1 

POSITION  ERRORS  DUE  TO  10  kt  AND  50  kt  WINDS 


10  kt 

UINO 

SO  kt  WIND 

TIM 

(•in  »9c) 

HORTH  POSITION 
ERROR  (na) 

EAST  POSITION 
ERROR  (na) 

NORTH  TOSITION 
ERROR  (na) 

East  position 

ERROR  (na) 

2  00 

1  0  ,  0  j 

1  0 

1 

0 

SO  * 

0.2*.  ;  0.021 

1 

0  098 

20:00 

0.0080  1  0.0010  i 

•0.0010  1 

-0  0020 

22:aO 

-0.20  j  0.21 

-1.6 

1  .0 
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